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ABSTRACT 

We investigate the relation between star formation rate (SFR) and gas surface densities in Galactic 
star forming regions using a sample of young stellar objects (YSOs) and massive dense clumps. Our YSO 
sample consists of objects located in 20 large molecular clouds from the Spitzer cores to disks (c2d) and 
Gould's Belt (GB) surveys. These data allow us to probe the regime of low-mass star formation essen- 
ce ■ tially invisible to tracers of high-mass star formation used to establish extragalactic SFR-gas relations. 

We estimate the gas surface density (E gas ) from extinction (Ay) maps and YSO SFR surface densities 
. (S§fr) from the number of YSOs, assuming a mean mass and lifetime. We also divide the clouds into 

evenly spaced contour levels of Ay, counting only Class I and Flat SED YSOs, which have not yet 
migrated from their birthplace. For a sample of massive star forming clumps, we derive SFRs from the 
total infrared luminosity and use HCN gas maps to estimate gas surface densities. We find that c2d and 
GB clouds lie above the extragalactic SFR-gas relations (e.g., Kennicutt-Schmidt Law) by factors up 
to 17. Cloud regions with high E gas lie above extragalactic relations up to a factor of 54 and overlap 
with high-mass star forming regions. We use 12 CO and 13 CO gas maps of the Perseus and Ophiuchus 
clouds from the COMPLETE survey to estimate gas surface densities and compare to measurements 

■ from Ay maps. We find that 13 CO, with the standard conversions to total gas, underestimates the Ay- 
Oh" based mass by factors of ~4-5. 12 CO may underestimate the total gas mass at E gas >200 Mq pc~ 2 by 
q , >30%; however, this small difference in mass estimates does not explain the large discrepancy between 
!— < ■ Galactic and extragalactic relations. We find evidence for a threshold of star formation (E t h) at 129±14 
"qq ' Mq pc~ 2 . At Sgas > S t h, the Galactic SFR-gas relation is linear. A possible reason for the difference 

between Galactic and extragalactic relations is that much of E gas is below E t h in extragalactic studies, 
which detect all the CO-emitting gas. If the Kennicutt-Schmidt relation (Esfr oc £ ga t) and a linear 
relation between dense gas and star formation is assumed, the fraction of dense star forming gas (/dense) 
increases as ~ E° a ^.. When E gas reaches ~300E t h, the fraction of dense gas is ~1, creating a maximal 
starburst. 

■ Subject headings: stars: formation — infrared: stars — ISM: dust — ISM: clouds — galaxies: ISM 

■ _ 

1. introduction E gas me asurements for the sample of galaxies in lKennicuttl 

o. 

, lar medium (ISM) control star formation is an important 
prerequisite to understanding galaxy evolution. A robust 



OO 



< 



TT , , , u • i • 1 1 • i i i (1199 8b) were fitted to a power law with N =1.4, which is 

Understanding now physical processes m the mterstel- , . ,, . _ . . ., , „ 

known as m the Kennicutt-Schmidt law ■ 



measurement of the relation between the star formation „ / c \ 14±015 



rate surface density (E SFR ) and the surface density of cold s sfr = (2.5±0.7) x 10 ( J (M Q yr kpc ) 



gas (E gas ) is of vital importance for input into theoretical 



(2) 



X 

H models of galaxy evolution. Since it ig Qnl &n assumption that there is only one 
The idea that there should be a relation between the rdaUm that regulates how gas is forming starSi we re _ 
density o f star fo rmaton and gas density was first pro- frain from caUing tHg & w and instead refer tQ it 
posed by i SchmidH Oa59|). Schmidt investigated this rela- ag & SFR _ gas re i ation5 or as the Kennicutt-Schmidt rela- 
tion now known as the Schmidt law , assuming that it tion w hen re f e rring to equation 2 specifically. Sever al au- 
shou d be m the form of a power law and suggested that thors (gario7][l9 92: Elmegreed [1991 Won g. fc Blitz! 12003 
the density of star forma tion^as proportional to gas den- iRrumholz fc argue that there is a simple expla- 
sity squared. |Kjpi£utj] U22SH) measured the global or nation for the Kennicutt-Schmidt relation: if the SFR is 
disk-averaged Schmidt law m a sample of spiral and star- proportiona i to the gas mass divided by the time it takes 
burst galaxies using the projected star formation and gas tQ CQnvert thg gag ^ starg &nd if we take ^ timescale 
surface densities (E gas ) m the form: to be the free _ faU tim6j then f K o.5 and ^ gFR a p i.5 



E S FR oc E^ s , (1) 



gas- 

Taking the scale height to be constant, p oc E, and this in 
J s as ' 1 1 ' turn gives the Kennicutt-Schmidt relation (to the extent 

where N is the power law index. The global SFR and that 1.4 ± 0.15 = 1.5). A variety of observational meth- 
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ods have been used to investigate this relation in different 
types of galaxies and on different scales. 

There have been many observational studies of SFR 
gas re lations on either glob al scales (IKennicutt 



1998b) or using either radial ([Martin fc Kennicuttl 
Wong fc BUtll2002tlBoissier et alJl2003tlHever et al 
U 



1989 



2001 



-■ ■20041: 

Schuster et all 120071) or point-by- 



Komugi et all l2005t 

point measurements ()Kuno et al.1 19951: IZhang et al.ll200lD 
that find values of N ranging from 1-2. Recently, there 
have been studies that measure star formation and gas 
content of galaxies on spatially resolved scales of ~0.2- 
2 kpc. These studies have fo und p ower law indices of 
N w 0.8 - 1.6 (IKennicutt et al.ll2007t pThilker et alj|2007 



Bigie l et al 



Verlev et al.l 



2008; IBraun et al.l 120091: IBlanc et al 



2009; 



2010). The study bv lBigiel et al.l ()2008D used 



a sample of 18 nearby galaxies to derive a spatially resolved 
relation on —750 pc scales. They found a linear relation 
between Ssfr and the molecular gas surface density over 
a range of 3-50 M pc~ 2 : 



J SFR. 



10 



-2.1±0.2 



10M Q pc" 



1.0±0.2 



(3) 

Other spatially resolved studies were based on measure- 
ments idoj3£jji_jJ 1 _shi^ pc: 
M51 (IKennicutt et alj |2007t [Blanc et "aTll200l NGC 7331 
([Thilker et alll2007l) . M31 (|Braun et al.ll2009D . and M33 
( Verlev et al.ll2010D . Since the global study of iKennicuttl 
|l998bj) and spatially resolved study of lBigiel et al.l (|2008| ) 
obtain results for large samples of galaxies, we use these 
studies as a baseline for comparison to our work. It is 
evident that sensitivity of N to systematic variations in 
methodology (e.g., data spatial resolution, SFR tracers, 
power law fitting method) accounts for some of the differ- 
ences in the derived star formation power law index, but 
the underlying physical reasons for the variations in the 
S FR-gas relations remain an open, challenging question. 

iKrumholz et a l. (2009) revisited the SFR-gas relation, 
considering the dependence on atomic and molecular 
components of £ gas , metallicity, and clumping of the 
gas. Their analysis produces a SF-gas relation that rises 
steeply at low £ gas , where the gas is mostly atomic, is 
nearly linear i n the regime where n o rmal spiral galax- 
ies are found (IKennicutt et al.l l2007t iBigiel et al.l l2008t 
Blan c et al.l 120091 ). and increases superlinearly above 85 
Mq pc~ 2 . Measurements made in these studies, how- 
ever, are limited to hundred parsec scales or more and are 
not directly comparable to the size of individual molecu- 
lar clouds or dense clumps where stars form. While these 
studies have all looked at the extragalactic SFR-gas rela- 
tion, there has been little work until recently investigating 
this relation locally in the Milky Way. 

Surveys of nearby molecular clouds in the Milky Way 
using Spitzer imaging have provided large statistical 
samples of young stellar object (YS O) candidat e s (e.g ., 
Allen et al. 2 010, i n preparation, lEvans et"aT1 ([20091) . 
iForbrich et all (I2009D . iRebull et all (12010D k These sur- 
veys have allowed us to directly count the number of low- 
mass stars that are forming and estimate SFRs. These 
data also allow us to trace the low-mass star formation 
regime essentially invisible to tracers, such as emission 
in Ha, ultraviolet, far-infrared (FIR), and singly ionized 



oxygen, used to establish extragalactic SFR-gas relations. 
Since these tracers only probe the rate at which massive 
stars form, a stellar Initial Mass Function (IMF), extrap- 
olating down to low stellar masses, must be assumed to 
obtain a SFR. Thus, these SFR estimates are very sensi- 
tive to the IMF slope and distribution on the low-mass 

end; 

lEvans et al.l ([2009D compared extragalactic observed 
SFR-gas relations to total molecular cloud measurements 
from the Spitzer c2d survey. They found that Galactic 
clouds lie abo ve the SFR-gas relations predicted b y extra- 
galactic work ([Bigiel et al.l 2008; K ennicuttlll998bl) and lie 
slightly above the ext rapolated relatio n from a study of 
massive dense clumps (|Wu et al.ll2005h : 



SFR - 1.2 x 10" 



M 



dense 



Mr; 



(M yr- 1 ). 



(4) 



This result suggests that studying SFR-gas relations in 
our Galaxy may be useful for interpreting star formation 
observed in nearby and high-z galaxies. On the high- 
mass end of the spectrum, a large survey of massive dense 
clumps bv lWu et al.1 ([2010), provides a sample that can be 
directly compared to extragalactic star forma tion tracers. 

In t his paper, we extend the comparison bv lEvans et al.l 
( 2009) by combining the 7 c2d clouds and 13 clouds from 
the GB survey. Regions of high-mass star formation 
from a su r vey of —50 massive dense Galactic clumps from 
IWu et all ([2010D provide an extension to high-mass star 
formation regions. The layout of this paper is organized 
as follows. We discuss low-mass star formation in the c2d 
and GB clouds and describe how £ gas is derived from ex- 
tinction maps and estimate SFR surface densities (Esfr) 
by YSO counts in Sections 12.11 and 12. 2[ respectively. In 
Section 12.31 we separate clouds into evenly spaced con- 
tour intervals of £ gas , measuring the Ssfr and S gas in 
these intervals. Section [3] discusses the differences between 
Galactic and extragalactic gas and SFR surface density re- 
lations. Ay and CO measurements of S gas are compared 
in Section 13.11 We investigate whether massive star form- 
ing regions behave differently from low-mass star forming 
regions in Section 13.21 The effects of averaging over whole 
galaxies (kpc scales), including both star forming gas and 
diffuse molecular gas, on the SFR-gas relation measured in 
extragalactic studies are discussed in Section 1331 Finally, 
we summarize our results in Section |4j 

2. LOW-MASS STAR FORMING REGIONS FROM Spitzer 
C2D AND GOULD'S BELT SURVEYS 

The cores to disks (c2d) Legacy project included 5 
large clouds: Serpens (Ser), Perseus (Per), C hamaeleon II 
(Cha II), Ophiuchus (Oph), and Lupus (Lup) ([Evans et al.l 
2009). Because the Lup 'cloud' is really composed of 
several separate clouds, we divide them in this study 
by name: Lup I, III, and IV, and obtain a total of 7 
clouds. The Gould's Belt (GB) Legacy project (Allen 
et al. 2010, in preparatio n) includes 13 larg e clouds: 
IC 5146E and IC 5146NW (lHarvev et al.ll2008l) . Corona 
Australis (CrA), Scorpius (Sco), Auriga (Aur), Auriga 
North (AurN), Serpens- Aquila ( Ser-Aqu), Musca (Mus), 
Cepheus (Cep) dKirk et al.ll2009h . Cha I and III, and Lup 
V and VI. These 20 clouds span a large range of masses, ar- 
eas, and number of YSOs (see Table[T|). The term "large" 
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was used in the c2d study to distinguish them from the 
sample of small clouds and cores th at were biased to - 
ward regions known to have dense gas ()Evans et a l. 2003) . 
The "large" clouds are thus suitable for statistical analyses, 
such as those presented here, but they are actually small 
compared to the Orion cloud or many clouds in the inner 
Galaxy. 

2.1. Estimating £ gas from Extinction Maps 

We derive cloud masses (M gas , cloud) and mean sur- 
face densities (£ gas , cloud) from the extinction maps, which 
were produced from a combination of 2MASS and Spitzer 
data, ranging from 1.25/im to 24/im. In this wave- 
length range, the spectral energy distributions (SEDs) of 
sources classified as stars provide measurements of the 
visual extinction (Av) along lines of sight through the 
clouds (Evans et al. 2007; Huard et al. 2010). Line-of- 
sight extinct ions were determined by fitt ing the SEDs, 
adopting the Wcingart ner fe Draind (|200ll) extinction law 
with Ry=Av I E(B — V)=5.5. Extinction maps were con- 
structed by convolving these line-of-sight measures with 
uniformly spaced Gaussian beams. The c2d team ob- 
served "off-cloud fields" for four of the molecular clouds: 
Chamaeleon, Perseus, Ophiuchus, and Lupus. The line- 
of-sight extinction measurements from these off-cloud 
fields suggested Ay calibration offsets of 1-2 magnitudes; 
therefore, in constructing the maps for these four clouds, 
Evans et al. (2007) subtracted these calibration offsets. 
Since no off-cloud field had been observed for Serpens, 
they used a weighted mean of the Ay calibration off- 
sets to correct the calibration of their Serpens extinction 
maps. No off-cloud fields were observed for the GB sur- 
vey. Further analysis of the fitting of the line-of-sight ex- 
tinctions demonstrates that the inferred calibration off- 
sets strongly depend on which wavebands had detections 
(Huard et al. 2010). For example, sources with only 
near-infrared (2MASS) detections may suggest no calibra- 
tion offset, while sources with only mid- infrared (Spitzer) 
detections show greater calibration offsets, perhaps as 
high as 2-3 magnitudes. Th is finding suggests that the 
Wcingar tner fc Draind (|200lh extinction law does not ac- 
curately characterize the reddening through the full range 
1.25-24 yitm spectral range. For if it did, the inferred Ay 
calibration offsets should be independent of the detected 
wavebands. For this reason, the extinction maps deliv- 
ered by the GB survey make use of the catalogued line- 
of-sight extinctions with no correction for potential cali- 
bration offsets, and, for consistency, they suggest that the 
previously adopted Ay calibration "offsets" of 1-2 mag- 
nitudes be added to the c2d extinction maps (Huard et 
al. 2010). After revising accordingly the extinctions in the 
clouds mapped by c2d, we find that the gas masses, and 
thus the cloud surface densities, are ~20-30% greater than 
those previously published by Evans et al. (2009). The ex- 
tinction maps used in this study probe to h igher Ay (up 
to 40 mag) than some previous stud ies (e.g.. iPineda et all 
( 20081) ; lLombardi et alj (|2008L l2010h ) due to the inclusion 
of both 2MASS and mid-IR Spitzer data. 

In order to compute the M gaS;C i ou( j, we chose extinc- 
tion maps with 270" beams for all clouds. We base this 
choice on the best resolution map available for Ophiuchus, 



which is limited in resolution due to a large extended 
region of high extinction with relatively few background 
stars detected. M gaS:C i ou( j and S gas , cloud were calculated 
by summing up extinction map measurements and con- 
verting to the column density using the relation NgJAv 
= (l mGC ntiV))- 1 = 1.37xl0 21 cm~ 2 mag" 1 (|Draind 
120031 ) for a lWeingartner fc Draind (|200ll ) R v =5.5 extinc- 
tion law, where C ext (V) — 6.715 xlO -22 cm 2 /H from the 
on-line tables^, using equation [5] and [51 respectively. The 
uncertainties in M gaS)C i ou( j are computed from maps of ex- 
tinction uncertainty, which account for the statistical pho- 
tometric uncertainties, but not systematic uncertainties in 
using the extinction law calibration. 

We compute M gaSjC i ou d by summing up all pixels 
(J2 Av) above Ay = 2 in all clouds except for Serpens 
and Ophiuchus which are covered by the c2d survey com- 
pletely down to Ay = 6 and 3, respectively. M gaSiC i ou( j is 
then: 



M„ 



cloud 



/im H (l.086C e: 
1.58 x 10" 36 x 



(v)Y 

Ay 

mag 



-pixel \ 



x A 



pixel 



cm 



(M X5) 



where the mean molecular weight (p) is 1.37, the total 
number of hydrogen atoms is N(H)=N(Hi) + 2N(H 2 ), and 
we assume a standard molecular cloud composition of 63% 
hydrogen, 36% helium, and 1% dust, mn is the mass of hy- 
drogen in grams, the area of a pixel in square cm (Api xe i) in 
the extinction map is (tt/180/3600) 2 D(cm) 2 R(") 2 , where 
R(") is the pixel size in arcseconds, and ^4 c ioud is the area 
of the cloud measured in square pc. We divide M gaS!C i ou d 
by the area to obtain £ gas , cloud for each cloud: 



cloud 



cloud 



-'gas, cloud 



pc^ 



15 



Ay 
mag 



(M pc" 2 )(6) 



(M pc" 2 ).(7) 



Measured cloud properties for c2d and GB clouds within 
a contour of Ay > 2 or Ay completeness limit are shown 
in Tabled] 

2.2. Estimating Star Formation Rates from YSO Counts 

We estimate the SFR from the total number of YSOs 
(A^YSO.tot) contained in an area where Ay > 2, as de- 
scribed in § 12.11 We assume a mean YSO mass ((Myso)) 
of 0.5±0.1 Mq, where the mean estimated error in mass 
is de rived from the mas s distribution of YSOs in Cha II 
from lSpezzi et all (|2008l ) . The mean YS O mass is consis - 
tent with IMF studies b vlChabrierl (|2003l) : iKroupal (|2002D : 
iNinkovic fc Traikovskal (|2006l) . We also assume a period 
for star formation (iciassii) of 2±1 Myr, based on the esti- 
mate of the elapsed time between form ation and the end of 
the Class II phase (jEvans et al.ll2009D . This SFR assumes 
that star formation has been continuous over a period 
greater than tciassii- All clouds have Class III objects, in- 
dicating that star formation has continued for longer than 
iciassii- The SFR measured in this way could be underes- 
timated or overestimated in any particular cloud, but over 
an ensemble of 20 clouds, it should be the most reliable 



SFR indicator available because no extrapolation from the 

5 Tables available at http://www.astro.princeton.edu/~draine/dust/dustmix.html 
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massive star tail of the IMF is needed. We base our er- 
ror estimates by choosing the largest error from either the 
systematic error, combined in quadrature from mean YSO 
mass and period of star formation, or the Poisson error 
from YSO number counts. 



Sspr = A 7 "-* 



YSCXtot 



/ (M- 



YSO/ 



V M d 



^cloud 

kpc 2 



, ^ClassII , 



(Moyr-^pc- 2 ) 



(8) 



Table[T]lists values for clouds within a contour of Ay > 2 
for all c2d and GB clouds. We show our estimated 
Sgas, cloud and Ssfr for the c2d and GB clouds in Figure[T] 
Sgas, cloud ranges from ~50-140 Mq pc -2 , and Esfr ranges 
from ^0.4-3.4 Mq kpc -2 yr -1 . We use these units for con- 
venience in comparing to the extragalactic relations. 

We compare the observations to the predicted values for 
Ssfr using £ gas , cloud that we calculate for the c2d and GB 
clouds. We plot these extragalactic relations in Figure Q] 
and will include them in all the following SFR-gas rela- 
tion figures. The solid lines represent the regime where 
they were fitted to data and the dashed lines are extrapo- 
lated relations spanning the range of E gas . The blue line 
is from disk-averaged or global SFR measurements based 
on Ha emission and the total (H i+CO) gas surface den- 
sities hi_a_sainple_ofjiormal spirals and starb urst galaxies 
from Kennicuttj (|l998bl) . The red line is from lBigiel et al.l 
(2008|), who made sub- kpc resolution measurements in a 
sample of spiral and dwarf galaxies using SFRs based on 
a combination of Spitzer 24/im and GALEX UV data and 
use CO measurements to obtain a relation for H 2 gas sur- 
face density. Both of these studies trace either obscured 
(24/im) or unobscured (Ha and UV) massive star forma- 
tion and are blind to regions of low-mass star formation 
that we are measuring in this work. We also compare to 
the theoretical total (H I + CO) gas and SFR relation of 
iKrumholz et al.l (|2009T> (orange solid line). This predic- 
tion takes into account three factors: the conversion of 
atomic to molecular gas, metallicity, and clumping of the 
gas. For our comparisons, we choose galactic solar metal- 
licity and a clumping factor of 1, which corresponds to 
clumping on 100 pc scales. We include data points for 
the Taurus molecu lar cloud, including YSO counts from 
iRebull et ail (I2010D . S gas from a 2MASS extinction map 
dPineda et all |2010|). and the tot al 13 CO and 12 CO gas 
mass from iGoldsmith et al.l (|2008l) . 

If we take the average £ gas , cloud defined as the total 
^gas, cloud divided by the total area (^4 c ioud) in pc 2 , we 
would find that the average molecular cloud in this study 
has a surface density of 91.5 M Q pc~ 2 and a Ssfr of 
1.2 M Q kpc~ 2 yr -1 . Taking this average £ gas , cloud and 
calculating what the extragalactic relations would predict 
for the average cloud SFR surface de nsity, we would ge t 
0.13, 0.07, a n d 0.03 M G k ^c'V" 1 forlKennicuttl (|1998bD . 
iBigiel et al.l (|2008l ). and IKrumholz et al.l (|2009D . respec- 
tively. The observed values exceed the observed extra- 
galactic Ssfr predictions by factors of ~9-17 and the the- 
oretical prediction by a factor of ~40. While the star for- 
mation surface density, Esfr of 1.2 Mq kpc -2 yr , seem 
high, the clouds fill only a small fraction of the local square 
kpc. From Table 1, the total SFR is 780.5 M Myr" 1 . 



If we remove the IC5146 clouds, which are more distant 
than 0.5 kpc, the SFR within 0.5 kpc is 748 M Q Myr _1 
or 7.5 x 10 -4 M yr . Extrapolated to the Galaxy with 
a star-forming radius of 10 kpc, this would amount to 0.3 
M Q yr -1 , less than the rate estimated for the entire galaxy 
of 0.68 to 1.45 M Q vr- ^Robitaille k Whitnevll2010D . This 
local, low-mass star formation mode thus could account for 
a substantial, but not dominant, amount of star formation 
in our Galaxy. 

2.3. Estimating S gas and Ssfr for the Youngest YSOs 
as a Function of Ay 

The last section gave us estimates over the whole molec- 
ular cloud including all YSOs in each c loud. Earl y 
work surveying large areas of clouds (e.g., lLadal (11992D ) 
suggested that star formation is concentrated in regions 
within molecular clouds in regions of high densities (n ~ 
10 4 cm -3 ). The c2d and GB studies of many whole 
clouds have clearly established that star formation is not 
spread uniformly over clouds, but is concentrated in re- 
gions at high exti nction. Furthermor e, the youngest YSOs 
and dense cores (|Enoch et al.ll2007D are the most highly 
concentrated at high Ay (Evans et al. 2009, Bressert et 
al. 2010, submitted). Older YSOs can leave their original 
formation region or even disperse the gas and dust. Tak- 
ing the average velocity dispersion of a core to be 1 km 
s -1 , a 2 Myr old YSO could travel ~2 pc, roughly the av- 
erage radius of a cloud in this study. We therefore apply a 
conservative approach and only estimate th e SFRs using 
the yo ungest Class I or Flat SED YSOs (see lGreene et al.l 
( 1994J) for the definition of classes) that have not yet mi- 
grated from their birth place. To classify YSOs as Class 
I or Flat SEP, we use the ex tinction corrected spectral 
index from lEvans et all (|2009l ) for the c2d clouds and the 
uncorrected spectral index for the GB clouds (Table [T|). 
These two classes of YSOs have timescales of 0.55± 0.28 
and 0.36± 0.18 Myr, respectively (Allen et al. 2010, in 
preparation) . 

In order to measure Ssfr and £ gas for the youngest 
YSOs, we divide the clouds into equally spaced contour 
levels of Ay or S gas ,con and measure the SFR, mass 
(M con ), and area (A con ) enclosed in that contour level. 
The contour intervals start from the extinction map com- 
pleteness limits (§ I2.ip and are spaced such that they are 
wider than our map beam size of 270" as shown in Fig- 
ure [2j We compute the gas surface density (E gas ,con) in 
the same way as in equation [5J but this time using only 
the mass (M g as,con) and area (A con ) enclosed in the Ay 
contour region. 



A. 



(Mq pC" 



(9) 



Mq J \ pC^ 

If there are no YSOs found in the contour region, we 
compute an upper limit to the SFR by assuming that there 
is one YSO in that region. The upper limits are denoted 
by the asterisks in Table [5J We estimate the uncertainties 
in the both the SFR and Ssfr by choosing the largest 
error: either the systematic or Poisson error from YSO 
counts. 

2.3.1. MISidentified YSOs from SED FITS (MISFITS) 
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The c2d and GB surveys have classified YSOs based on 
the SED slope from a fit to photometry between 2{im and 
24/im (Evans et al. 2009; Allen et al. 2010, in prepara- 
tion). However, we find that some of the Class I and Flat 
SED YSOs are not clustered and lie farther from the ex- 
tinction peaks than expected for their age. Most of these 
suspicious objects are found, on average, to lie at Ay —6 
magnitude. If these Class I and Flat SED objects are true 
young YSOs, they are more likely to be cent rally concen - 
trated toward the densest regions in a cloud (|Ladalll992h . 
Class I and Flat SED YSOs should be associated with a 
dense, centrally concentrated, envelope of gas. We there- 
fore performed a follow up survey of these suspicious ob- 
jects for a subset of c2d and GB clouds using the Caltech 
Submillimeter Observatory (C SO). Our work wa s mot i- 
vated by a study performed by Ivan Kempen et alj (|2009l ) . 
who mapped HCO + J— 4-3 using the James Clerk Maxwell 
Telescope and found that 6 previously classified Class I 
YSOs in Ophiuchus had no detections down to 0.1 K. 

With high effective an d critical densities, n~10 4 and 
-10 6 cm" 3 (|Evansl |1999( ). the HCO+J=3-2 transition 
provides a good tracer of the dense gas contained in pro- 
tostellar envelopes. We observed Class I and Flat SED 
YSOs at the CSO from the Aur, Cep, IC5146, Lup, Oph, 
Per, Sco, Ser, and Ser-Aqu molecular clouds using the 
HCO+ J=3-2 (267.557620 GHz) line transition as a dense 
stage I gas tracer to test if they are truly embedded YSOs. 
A survey of all observable c2d and GB clouds with detailed 
results will be published in a later paper. 

Observations were made during June and December of 
2009 and July 2010 with an atmospheric optical depth 
(T225) ranging from 0.05-0.2. We observed each source us- 
ing position switching for an average of 120 seconds on and 
off source. If a source was detected, we integrated until we 
reach a signal-to-noise of > 2a (most sources have > 3a 
detections). Using an average main beam efficiency (77) of 
0.80 and 0.61 for December 2009 and July 2010 observing 
runs, respectively, we compute the main beam tempera- 
tures and integrated intensities of detected sources. The 
results are shown in Table [3] 

For this paper, we observed a total of 98 suspicious 
sources, 45 Flat SED and 53 Class I sources. We find that 
74% (73/98) of the observed sources are not detected in 
HCO+ J=3-2. Out of the 42 Flat SED sources, we detect 
only 3, but we detect 42% (22/53) of the Class I sources. 
The YSO MISFITS are a small fraction of the total num- 
ber (3146) of YSOs or Class I plus Flat sources (681) in 
the c2d and GB studies, but they could bias the statistics 
upward at low gas surface densities. The undetected MIS- 
FITS may be background galaxies or later stage YSOs, and 
we will explore this in more detail in a later paper. Fig- 
ure [H shows the distribution of Class I (red filled circles), 
Flat SED YSOs (yellow filled circles) and non detected 
MISFITS, indicated by the open stars on the Per cloud 
Ay map. MISFITS that we do not detect in HCO+ J=3- 
2 are removed from the sample when we measure S gaSjCon 
and Ssfr in Ay contours (Section [ 



cloud or separate cloud component (see Section \2§ . In 
Figure [3l we show the S gaSjCon and Ssfr densities for 
both Class I and Flat SED sources (green and magenta 
stars) and upper limits for each class (green and magenta 
inverted triangles) that we measured in contour regions 
described in § 12.31 with extragalactic observational rela- 
tions over-plotted. A wider range in both S gas (~45- 
560 M Q pc" 2 ), and S SF r (-0.03-95 M Q kpc~ 2 yr" 1 ) 
are found for contour regions compared to the total cloud 
measurements. W e note that the points for Sco and Cep 
(Kir k et al.l 12009) clouds are obtained by co-adding the 
separate cloud regions using the same contour intervals. 
Since these points sample regions with non-uniform Ay, 
they only provide an estimate of Ssfr and S gas . These 
points lie at Ssfr < 6 M kpc -2 yr^ 1 and high S gas >330 
M Q pc- 2 . 

We compare our YSO contour results to extragalactic 
relations and find that most points lie well above the ex- 
tragalactic relations. Excluding upper limits, the mean 
values of Ssfr and S gaS)COn are of 9.7 M kpc~ 2 yr _1 and 
225 M pc~ 2 , respectively. Evaluated at this mean gas 
surface density, the extragalactic relations under-predict 
Ssfr by f actors of —21-54 . The mean YSO contour lies 
above the iKrumholz et al.1 (|2009| ) extragalactic SFR-gas 
relation prediction by ~2 orders of magnitude. We ex- 
plore the differences between the Galactic and extragalac- 
tic SFR-gas relations in Section [3] 

3. WHY ARE GALACTIC SFR-GAS RELATIONS DIFFERENT 
FROM EXTRAGALACTIC RELATIONS? 

The differences between our findings on Galactic scales 
and the extragal actic relations, bo th on global or disk- 
averaged scal es dKennicuttl Il998bl) and scales of hun 
Kenrricutt et all 120071 iThilker et al " 
iBlanc et al.1 120091 : iBraun et al. 



dreds of pc (H 
Bigi el et all l"2008; 



2007 



2009; 



2.4. Results: The Youngest YSOs as a Function of S gas 

After removing the MISFITS from our Class I and Flat 
SED YSO sample, we show the number, Mg aSjCon , S gas con , 
SFRs, and Ssfr in Table [2] for all contour levels in each 



Verlev et al.l |201Q[) , might be explained in the following 
ways. Firstly, using 12 CO to measure the H2 in galax- 
ies might give systematically different S gas than do Ay 
measurements (Section |3.1[) . Secondly, the local c2d and 
GB clouds are forming low-mass stars; since extragalactic 
SFR tracers respond only to massive stars, the two star 
forming regimes might behave differently. In Section 13.21 
we will investigate whether massive star forming regions 
agree with the extragalactic SFR-gas relations and if they 
vary from low-mass star forming regions. Finally, aver- 
aging over whole galaxies on scales of hundreds of pc in- 
cludes both gas contained in the parts of molecular clouds 
that are forming stars and diffuse molecular gas that is 
not forming stars (Section 13.31) . A local example of this 
is a st udy of the Taurus molecular cloud; iGoldsmith et al.l 
(2008), found a large amount of diffuse 12 CO at lower gas 
densities where no young stars are forming. Extragalactic 
studies averaging over hundreds of pc scales would include 
this gas, causing an increase in the amount of CO flux that 
is being counted as star forming gas. 

3.1. The use of CO versus Ay to determine S gas 

Since extinction maps are direct probes of S gas , they 
provide the best measure of the total gas and are opti- 
mal for use in determining the S gas of molecular clouds. 
However, Ay maps are not easily obtainable in extragalac- 
tic studies, which instead employ CO maps, particularly 
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12 C0 J=l-0, to determine E gas of molecular hydrogen. 
Since the molecular hydrogen (H2) rotational transitions 
require high temperatures not found in the bulk of molec- 
ular clouds, other tracers of dense gas are used to estimate 
the amount of H2 . The next most abundant molecule with 
easily observable excitation properties in a molecular cloud 
is 12 CO J—1-0. In this study, we want to explore how 
well CO traces Ay as a function of M gas or S gas . We can 
directly test this in two galactic clouds, Perseus and Ophi- 
uchus, which both have 12 CO J=l-0 and 13 CO J=l-0 
maps from the Five College Radio Astronomy Observatory 
(FCRAO) Coordinated Molecular Probe Line Extinction 
Thermal Emission (COMP LETE) Survey of Star Forming 
Regions (jRidge et al.ll2006l) . 

In order to directly compare the CO maps from the 
COMPLETE survey to the Ay maps in this study, we 
interpolate the CO data onto the Ay map grid with a 
pixel size of 45". We integrate the publicly available CO 
data cubes over the velocity range from 0-15 km s _1 
to create moment zero maps of integrated intensity de- 
fined as: Ico( x ,y) = J T m b(x,y, z)dV K km s _1 , where 
T m b = T^/r] m b is the main beam brightness temperature 
defined as the antenna temperature (Tj£) divided by the 
main beam efficiency (r) m b) of 0.45 and .49 for 12 CO 
J=l- and 13 CO J=l-0, respectively, from iPineda et all 
(2008]). Corresponding rms noise maps (cT( x ,y)) were con- 
structed by calculating the standard deviation of intensity 
values within each spectroscopic channel where no signal 
is detected. In order to determine the gas surface density 
of H2 (Eh 2 ), we must first compute the column density 
of H2. The column density of H2 is estimated from the 
12 CO map by using a CO-to-H 2 conversion factor (Xco), 
which is defined as the ratio of H2 column density to the 
integrated intensity (Aco = Nn^/Ico)- Similarly for the 
13 CO map, the column density is derived from the 12 CO 
and 13 CO maps, assuming LTE and an abundance ratio 
of H 2 -to- 13 CO. To compare S H . 2 from 12 CO and 13 CO to 
the S gas from Ay, we use only regions in the CO maps 
that have emission lines with positive integrated intensi- 
ties and line peaks that are greater than 5 times the rms 
noise. Our masses from extinction measurements used for 
this comparison are therefore slightly lower than those in 
Tables Q] and [2] by ~5% (See Tables |1 and [5]) . 

The Xco factor for 12 CO has been derived us- 
ing a variety of methods such as gamma ray emis- 
sion cause d by the col li sion of cosmic rays with 
hydr ogen dBloemen et all 1 19861). virial m ass meth- 
ods (jSolomon et all Il987t [Blitz et al] l2007t) . maps of 
dust emission from I RAS and assuming a con- 
stant dust-to-gas ratio (j Frerking et al .1 1982 ). e xtinc- 
tion maps from optica l star c ounts dDuvert et all [l986; 



iBachiller & Cernicharol 119861: [Linger et alj Il989h an d 
2MASS data (jLombardi et all [20061: IPineda et al.M2008D . 
and theoretically by the assumption that giant molecu- 
lar clouds are in gravita tional equilibrium ()Dickman et al.l 
119861: iHever et al.l l2001). All these studies find a range of 



Ago of 0.9-4.8xl0 2U cm" 2 K _i km" 1 s, but they were 
almost all restricted to regions with Ay < 6 mag. Stud- 
ies of extragalactic SFR-gas relations chose values close to 
the average galactic Aco measurements in t he literature: 
2.0xl0 20 cm 



2 FT 1 km " 1 dBigiel et al.|2008D) o r2.8xl0 20 
- 2 K- 1 km- 1 s from lBloemen et all (|1986l ) (|Kennicuttl 



ll998btlKennicutt et al.ll2007tlBlanc et al.ll2009l ). Since the 
goal of this study is to compare to extragalactic measure- 
ments, we choose a Ag o of 2.8 ±0.7xl0 20 cm~ 2 K _1 km -1 



s from iBloemen et al 



study of lKennicuttn fT998b). 



to be consistent with the 



We compute the column density of H2 from 12 CO mea- 
surements using the equation: 

Nh 2 (12 C o) = Xi2 CO x 7i2 CO (cm" 2 ). (10) 
This can be rewritten in terms of gas surface density: 



s h 2 ( 12 co) — (2m H x A H2 (i2 CO ) x A„\ 



1Q -33 



x toh x 



H 2 (i2CO) 



pixel) /A c \ QU d 
^pixcl 



A. 



loud 



pc^ 



(M Q pc" 



where we take the total number of hydrogen atoms 
is N(H)=2N(H2). We use this factor of two instead of 
the mean molecular weight of H2 (/in 2 =2.8, derived from 
cosmic abundanc es of 71% hydrogen, 27 % helium, and 
2% metals, e.g., iKauffmann et al.l (I2008D) to be consis - 
tent with the ex t ragala ctic studies of iKennicuttj (|1998bf) 
and iBigiel et al.l (|2008l ). This factor of 2 does not ac- 
count for helium, which is an additional factor of ^1.36 
(|Hildebrandlll983h . The errors in our gas surface density 
measurements include both the error from the rms inten- 
sity maps and the error in t he CO-to-H2 conversion factor 
from lBloemen et al.l (|l986l ). 

The top two panels of Figure [4] show the 12 CO inte- 
grated intensity versus the visual extinction derived from 
the 2MASS and Spitzer data for both Per (left) and Oph 
(right). We over-plot the conversion factor der ived from 
the gamma ray study of Bl oemen et al.l (<1986l) (dashed 
line). 12 CO is seen to correlate with Ay out to Ay ~7-10, 
where 12 CO starts to saturate and the distribution flattens 
out to higher Ay. A large difference in the amount of 12 CO 
integrated intensity produced relative to that predicted 
by Ai2 CO between the Per and Oph molecular clouds is 
seen. This may be due to higher opacity at higher Ay in 
Oph relative to Per. Extinction values around 10-20 mag 
are found to be essentially invisible to 12 CO. This figure 
demonstrates the non-linear, non-monotonic behavior of 
CO emission with Ay. 

12 CO, however, is not the most reliable tracer of star 
forming gas because of high opacity and varying 12 CO-to- 
H2 abundance due to photodissociation o r depletion on to 
dust grains. Studi es of molecular clouds (ICarpenter et al.l 
119951; IHever et al.lll996t iGoldsmith et al.ll2008D show that 
12 CO contains a significant diffuse component in the low 
column density regime Ay <4. 13 CO emission is a more 
reliable tracer of dense gas ranging from 1000-7000 cm~ 3 
than 12 CO because it is optically thin for most con- 
ditions within a cloud and 13 CO abundance variations 
are sma ll for densities <5000 cm" 3 and temperatures of 
>15 K (IBachiller fc Cernichardll986l: iDuvert et alJ ll986t 
IHever fc Laddlll995l: ICaselli et al.lll999D . 

To estimate the column densities from the 13 CO J=l-0 
integrated intensity maps we assume local thermodynamic 
equilibrium (LTE), optically thin 13 CO J=l-0, and that 
12 CO J=l-0 and 13 CO J=l-0 have equivalent excitation 
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temperatures. In order to derive column densities, we also 
determine an optical depth (ti3 CO ) and excitation temper- 
ature (T ex ) from comparison to the 12 CO J=l-0 line. We 
can derive this by assuming the 12 CO J=l-0 is optically 
thick: as t — s> oo, 



or in terms of surface densities: 



5.5 



In 1 



5.5 



r peafc,12C0+ 



(K), 



(11) 



where Tp ea k . 12010 is the 12 CO J =1-0 peak main beam 
brightness temperature which we measure on a pixel-by- 
pixel basis. 

The 13 CO J =1-0 optical depth is then 



3 co 



-In 1 



T, 



pcak : 13 CO 

5.3 



1 



exp 



(5.3/T CJ 



-0.16 



(12) 

where T pca k ; i3 C q is the 13 CO J =1-0 peak main beam 
brightness temperature measured in each pixel. We can 
then use the definition of 13 CO optic al depth and column 
density from lRohlfs fe Wilson! (|l996l ) to estimate the col- 
umn density of 13 CO: 



iVi3 CO = 2.6 x 10 1 



T13CO 



1 — exp 



3 CO 



(13) 



(cm 



,1 -exp(~ 5 - 3 / T - 

Certain regions near Ay peaks in the clouds are optically 
thick and are affected by 12 CO self absorption. In these 
regions, we cannot accurately determine T ex and therefore 
Ti3co- This problem affects ~1Q% of the pixels in Oph 
(Ay > 30 mag) and ~5% of the pixels in Per (A v > 20 
mag). Since most of the mass lies at low Ay, we mask 
out the pixels that have 13 CO self absorption and do not 
include them in determining the 13 CO column density and 
mass estimate discussed below. 

For comparison with extragalactic work, where clouds 
are not resolved, we derive 13 CO column densities by us- 
ing average spectra to determine T ex and ri3co using a 
T pGak i2 CO of 3 and 5.3 K and a T peak i3co of 0.7 and 1.5 K 
for Per and Oph, respectively. Comparing the two meth- 
ods, we find that using peak temperatures from average 
spectra or a constant T cx and Ti3co wm result in higher 
7Vi3co by a factor of ^2 at Ay < 10 over the pixel-by- 
pixel measurements. 

In order to convert the Ni3qq column density into H2 
column density, we use a H2-to- 13 CO abundance ratio. 
The H2/ 13 C O abundance r atio fo r the Per cloud was de- 
termined by IPineda et al.1 (|2008l ). who found a value of 
3.98±0.07xl0 5 for Ay < 5 mag. Dividing the cloud into 
separate regions, they found an average abundance ratio 
of 3.8xl0 5 . Other values found in the literature range 
from 3.5-6.7xl0 5 , with an average value of ~ 4xl0 5 us- 
ing bot h extinction maps (IPineda et al.l 20081) and star 
counts (iBachiller fc Cernicharolll986t iDuvert et al.lll986t 
lLanger et al.lll989l ). We adopt the average H2/ 13 CO ratio 
from the literature of (4±0.4)xl0 5 to convert 13 CO to H 2 
column densities using the relation: 



H 2 (i3CO) 



(4 ±0.4) x 10 5 iVi3 CO (cm" 2 ) 



(14) 



Sh 2 («co) = (2m H x N H2( i 3 



loud 



10 



-33 



x toh x 



H 2 ( 13 CO) 



pc z 



A, 



(M pc" 



where we choose a factor of two instead of the mean 
molecular weight in order to consistently compare to 12 CO. 
We show the 13 CO integrated intensity versus the visual 
extinction derived from the 2MASS and Spitzer data for 
both Per (left) and Oph (right) in the bottom two panels 
of Figure 2J The average i?2/ 13 CO ratio is shown by the 
dashed line. A turnover in 13 CO is seen at Ay ~7 (Per) 
and ~10 (Oph) that is likely due to an increase in optical 
depth. The amount of 13 CO in Per follows the average 
abundance ratio out to Ay ~5, but in Oph, the 13 CO 
integrated intensity is underproduced. 

To test how well CO traces Ay as a function of M gas or 
Sgas, we measure E gas densities in our Ay maps and £h 2 
in the CO maps in the overlapping area where there is a 
positive CO integrated intensity over 5 times the rms noise. 
In Figure 03 we plot the ratio of £h 2 an d S gas from Ay, 
which are effectively mass ratios, since the area measured 
is the same. The cyan squares and circles are points for the 
c2d and GB clouds (£h 2 ( 12 CO,c1ou<1), £H 2 ( 13 CO,cloud)) and 
the filled green and yellow squares and circles are mea- 
surements in contours of Ay using the same method as in 
SectionO(S H2 (i2co,con), S H2 (i3 C o,con)) for Oph and Per, 
respectively (Tables 2] and [5]) . A measurement for the 
Taur us cloud using both 12 CO and 13 CO above Ay = 2 
from ([Goldsmith et al.|[2008l ) is also shown (cyan triangle) . 
If CO traces the mass we find using extinction maps, we 
would expect the ratio of CO/ Ay mass to be of order unity 
as shown by the solid black line in Figure [3] For 12 CO, 
we find the total cloud measurement for Per to have £h 2 
of ~ 1.6S gas at Sgas^ 100 Mq pc~ 2 , but the ratio is close 
to unity within the errors. We find that 12 CO traces Ay 
relatively well in the Oph cloud out to ^200 M© pc~ 2 . At 
Sgas ^ 200M© pc -2 , 12 CO underestimates the Ay mass 
in both Per and Oph by ~ 30%, on av erage. 

Since 13 CO should trace de nser gas (jDuvert et al.| [l986: 
Bac hiller fc Cernicharol |1986[ ). we also explore how it 
traces Ay as a function of £ gas - We plot this on Fig- 
ure [5] for the total clouds (cyan points) and contour mea- 
surements (green points). We find a constant value of 
13 CO versus the surface density of extinction, but find 
that it underestimates £ gas by a factor of ~4-5 and lies 
below measurements of 12 CO by a factor of ~5, on average. 
The difference we find between 13 CO and H2 measured by 
Ay , could be due to the LTE method we used to compute 
13 CO masses or the assumption t hat there is a cons tant 
abundance of CO relative to H2. iHever et al.l ((2009) ex- 
plored the properties of galactic molecular clouds using 
13 CO emission and found that the assumption of equiva- 
lent excitation temperatures for both 12 CO and 13 CO in 
the LTE method may underestimate the true column den- 
sity of 13 CO in subthermally excited regions. As the 13 CO 
density increases, the J =1-0 transition becomes thermal- 
ized, and t he column d e nsity esti mates are more ac c urate . 
Also, both IHever et al.l (<2009D and lGoldsmith et al.l |2008) 
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found that if 13 CO-to-H2 abundance variations in LTE- 
derived cloud masses are not considered, they would un- 
derestimate the true column densities by factors of 2-3. 
Since we only include 13 CO emission greater than 5 times 
the rms noise, we are likely measuring gas that is thermal- 
ized with little abundance variation. Assuming a constant 
abundance ratio will therefore not account for the differ- 
ence between Ay and 13 CO masses. 13 CO, might therefore 
be a more consistent tracer of E gas , but it may underesti- 
mate the mass by factors of ~4-5, which can be corrected 
for. 

Variations in the CO-to-H2 conversion factor may im- 
pact the slope of the SFR-gas relations as measured by 
extragalactic studies. Since we are resolving molecular 
clouds, we cannot place constraints on gas de nsities lower 
than ~50 M© pc -2 , typical of spiral galaxies (|Bigiel et al.l 
2008). However if we consider the effects of using CO as a 
tracer of the total gas density, it underestimates the mass 
measured from A v by >30% for £ gas > 2OOM pc" 2 . This 
would effectively shift the extragalactic observed points to 
the right above 200 M© pc~ 2 . This sh i ft wou ld flatten 
the slope slightly in the fitted IKennicuttl (|1998bl) relation. 
These small factors, however do not explain the large dis- 
crepancy between the extragalactic relations and the much 
higher SFR in the local clouds, seen both in the whole 
molecular clouds and looking at the youngest YSOs as a 
function of £ gas . 

3.2. Do High-mass and Low-mass Star Forming Regions 
Behave Differently? 

By studying nearby molecular clouds, we can obtain the 
most accurate measurement of £ gas and £sfr> but it is 
regions of massive star formation that form the basis for 
extragalactic studies. Massive star forming regions are the 
only readily visible regions forming stars at large distances 
and thus are the only probes of star formation in distant 
regions in the Milky Way and external galaxies. We can 
measure Xsfr and S gas in individual massive star forming 
regions to see if there is better agreement with extragalac- 
tic SFR-gas relations. 

To investigate where individual regions of massive star 
formation fall on the SFR-gas relation, we use data from 
the molecular line survey of dense gas tracers i n > 5 
massive dens e ((n) ~ 10 6 cm " 3 , e.g. iPlume et al.| (|1997ft ) 
clumps from lWu et all (|2010l ). The IWu et alj (|2010D sur- 
vey measured clump sizes, virial masses, and dense gas 
surface densities using HCN J=l-0 as a tracer (Ehcn) at 
FWHM of the peak intensity. These are the sites of for- 
mation of clusters and massive stars. The most popular 
tracers of massive star formation include the ultraviolet, 
Ha, FIR, and singly ionized oxygen; however, due to high 
extinction toward and in these regions, we can use only the 
total IR luminosity to measure the SFR in these massive 
clumps. 

Since HCN J=l-0 gas has been shown to be tightly cor- 
related with the total IR lumin osity in clumps as long as 
Lir > 10 4 - 5 L (IWu et al.ll2005l).aswell as in both normal 
spiral and star burst galaxies ( Gao fc Solomon! l2004al lbT). 
we can use it to compare gas and star formation from the 
total IR luminosity in both the Milky Way and external 
galaxies. 

Shcn is calculated using the mass contained within the 



FWHM siz e (-Rhcn) of the HCN gas, following the meth- 
ods used in lShirlev et~ail (|2003l) : 



^HCN = 7T 1 X 



Mvi 
Mr. 



V PC 



(M© pc" 2 ), (15) 



where M v ; r is the virial mass enclosed in the source size 
at FWHM intensity. Uncertainties in the Ehcn are com- 
puted by adding in quadrature t he errors in the FWHM 
size and the mass as discussed in lWu et al.l (j2010D . 

We compute the SFR for massive dense clumps follow- 
ing extragalactic methods using the total infrared (IR) 
luminosity (Lir; 8-1000/^m) derived from the 4 IRAS 
bands. We assume the conversion SFRtr (M© yr _1 ) ~ 
2xlO- 10 L m (Lq) from IKennicuttl (jl998al) . S S fr ir are 
computed using the FWHM source sizes: 



J SFRi R 



SFR 



IB 



M; 



©yr" 



-ikpc- 2 ). 

n<ii 

The uncertainties in the EgFR IH density only include the 
error in FWHM size and a 30% error in SFR calibrations 
using the IR. In fact, the uncertainties are larger. The 
SFR calibration assumes that the observed far-IR emis- 
sion is re-radiate d by dust heated by O, B, and A stars 
(|Kennicuttl 1199 8a) . For low SFR, heating by older stars of 
dust unrelated to star formation can contaminate the SFR 
signal, causing an overestimate of the SFR. This is not a 
problem for the regions of massive star formation in our 
Galaxy, where the heating is certainly due to the recent 
star formation. A bigger issue for the Galactic regions is 
that the full Lir seen in the extragalactic studies is not 
reached unless the individual region forms enough stars 
to fully sample the IMF and is old enough that the stars 
have reached their full luminosity. These conditions may 
not be met during the time span of an individual massiv e 
clump dKrumholz fc Thompson] l2007t lUrban et al.|[2010h . 
In particular. lUrban et al.l (|2010t ) have calculated the ratio 
of luminosity to SFR in a simulation of a cluster forming 
clump. They find that L/SFR increases rapidly with time, 
but lies a factor of 3-10 below the relation in IKennicuttl 
(1998af) when their simulations end at times of 0.7 to 1.4 
Myr. Therefore, we may expect both large variations and a 
tendency to underestimate the SFR in individual regions. 
Unfortunately, despite these issues, Lir remains the best 
measure of SFR available to us at present in these regions. 
While it would be suspect to apply a correction factor 
based on the Urban et al. simulation, increasing the SFR 
of the regions of massive star formation by 0.5-1 order of 
magnitude would bring them into better agreement with 
the highest surface density points from the nearby clouds. 

For the sample of ~50 massive dense clumps, 42 sources 
have corresponding IR measurements. The resulting gas 
surface densities and Ssfr for the sample of massive dense 
clumps are shown in Figure [5] and Table [5] The relation 
between SFR and dense gas mass, Md en se(H2), for galac- 
tic clumps can be derived from (Lir/Lhcn(i-o)) = 911 ± 
227(K km s- 1 pc 2 )" 1 , and (M dense (g 2 )/L H CNn -o)) = 
11±2M (K km s- 1 pc 2 )" 1 from IWu et aH ((20051) . These 
relations can be combined with the IR SFR conversion 
from above to obtain a relation for £sfr ir and Ercn: 
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^SFRi 



1.66 x 10" 



^HCN(l-O) 

1M S pc" 2 



(Moyr-^pc" 2 ). 



Thi s equation is equivalent to that shown in IWu et al.l 
(|2005l ). which is the fit to both massive clumps and galax- 
ies. 

IWu et all (|2005f) found a decline in the linear L — 
-^hcn(i-o) correlation at Lir < 10 4 ' 5 Lq, where the clump 
is not massive or old enough to sample the IMF. Since the 
majority of the points with Lir, < 1O 45 L in Figure [5] 
lie off this relation, we include only massive dense clumps 
withLiR > 10 45 L Q . The resulting number of HCNJ= 1-0 
sources is 25 (Tabic [5]). The HCN J= 1-0 clumps are found 
to range from ~ 10 2 -4.5 x 10 3 M Q pc~ 2 in Ehcn and 
from 2-130 Mq kpc~ 2 yr _1 in £sFR m - These £ gas and 
SsfRir values are s imilar to t h ose of circumnuclear star- 
burst galaxies from IKennicuttl (jl998bf ). which range from 
~ 10 2 -6xl0 4 M Q pc" 2 and 0.1-9.5 x 10 2 M kpc" 2 yr" 1 
(see Figure [H]). The average HCN J=l-0 clump has 
Shcn(i-o) of (l-3±0.2) x 10 3 Mq pc" 2 and E SFRlR of 28±6 
Mq kpc~ 2 yr _1 . 

We also compare the relation we find for the massive 
dense clumps to known extragalactic relations in Figure [7] 
In this figu re we compare to th e H2 gas surface density re- 
lation from lBigiel et al.l (|2008|) and the total (H 1+H2) gas 
surface density relation from IKennicuttl (|l998bl ). We find 
that mo st of the points for the HCN J= 1-0 line lie above 
both the lBigiel et ail (<2008D and IKennicuttl (<1998bD extra- 
galactic relations wit h the aver a ge clum p ly ing a factor 
of ~ 5-20 above the IKennicuttl (jl998bl ) and iBigiel et all 
(2008) relations, respectively. 

In Figure [5] we plot the ratio of £sFR/£gas versus £ gas 
for c2d and GB clouds, YSOs, and massive clumps. We 
find a steep decline in Esfr and SsFR/E gas at around 
~100-200 Mq pc -2 . We identify this steep change in Ssfr 
over ^100-200 M Q pc -2 as a star forming threshold (Eth) 
between regions actively forming stars and those that are 
forming few or no low-mass stars. 

In Figure[9]we show points for the massive dense clumps, 
c2d and GB cl ouds, the youngest YSOs, and both the 
IWu et all (|2005h and extragalactic relations. We also show 
the range of gas surface densities for spiral a nd circum- 
nuclear starburst galaxies from the sample of IKennicuttl 



(fl998bh . E 



gas, con 



for Class I and Flat SED YSOs lie in- 
termediate between the regions wh ere spiral galaxies and 
starburst galaxies are found on the IKennicuttl (|1998bD re- 
lation. At S gas > S t h, the youngest Class I and Flat 
SED YSOs overlap with the massive clumps (Figure [5]) ■ 
Therefore, high-mass and low-mass star forming regions 
behave similarly in the EsFR _ E gas plane. The difference 
between extragalactic relations and c2d and GB clouds is 
not caused by the lack of massive stars in the local clouds. 
Also, the overlap with the massive clumps in Figure [^sug- 
gests that Lir provides a reasonable SFR indicator, as long 
as it exceeds 10 4 5 L Q , though an upward correction would 
produce better agreement. 

A steep increase and possible leveling off in Esfr at 
a threshold E t h ^100-200 M Q pc~ 2 is seen in both Fig- 
ures [5] and [51 We further constrain this steep increase and 
the possibility of Esfr. flattening at E gas > Eth, by ap- 
proximating it as broken power law with a steep rise that 



levels off in Section 13.2.11 



3.2.1. Star Formation Threshold 

In order to determine a robust estimate of £ t h, we n t the 
data using two models: a single power law (y = Nx + A, 
where y =logio Ssfr! x= logio E gas and a broken power 
law (z/i = N-lx + Ai;y 2 = N 2 x + A 2 ). We first fit Class I 
and Flat SED YSO points (£sFR,conj E gaSjCon ) and massive 
clumps (EsfRir , Shcn(i-o)) to a single power law. We do 
not include upper limits for YSOs or points for Sco and 
Cep clouds, which are co-added separate cloud regions and 
only provide a rough estimate of Esfr and E gas . The sin- 
gle power law fit yields N = 1.57±0.09 and A = -3.0±0.2, 
and a reduced chi-square (% 2 ) of 3.7 (84 dof). We fit the 
data for both YSOs and massive clumps to a broken power 
law for the range of S gas =50-250 M Q pc 2 . We minimize 
the total x 2 for the two segments of the broken power law 
using a simplex routine, which yields best fit parameters: 
Ni = 4.58 ± 0.5, A x = -9.18 ± 0.9, N 2 = 1.12 ± 0.07, 
A 2 = -1.89 ± 0.2 with a xl of 3.04 (82 dof). We at- 
tribute the large Xr to the scatter and large errors in the 
data, but since the Xr is ~18% lower for the broken power 
law compared to the single power law, we take it to be 
the best fit model. Equating the broken power law fits 
[yi — 1/2)1 we obtain a power law break at E t h= 129.2 
Mq pc 2 (Ay = 8.6) with a statistical 1-sigma deviation 
in x 2 °f = ' = 14 Mq pc 2 giving a range in E t h of ~115- 
143M Q pc 2 . Figure [TO] shows the broken power law fit 
(cyan and magenta lines), S t h, and the 1-sigma statistical 
range of S t h (dashed black vertical line and grey shaded 
region, respectively). The slope of the broken power law 
changes from a steep relation at S gas < S t h (slope of ^4.6) 
to linear relation (slope of ~1.1) at E gas > E t h- We note 
however, that variations in cloud distances will change this 
threshold slightly. One e xample is for Ser, w hich has a re- 
cent distance estima te by Dzib et al.l (I2010D of 415 pc and 
another estimate by IStraizvs et alj ( 19961 ) of 260 pc (used 
in this paper). If we use the larger distance of 415 pc, 
this would change our star formation threshold slightly to 
126±12 M© pc 2 (Ay = 8.4). 

This star forming threshold we find is in agreement 
with the threshold found in a study of local molecular 
clouds by Lada, Lombardi, & Alves (2010, subm itted) 
at A v - 7 or 116 M Q pc 2 . lEnoch et all (|2007l) also 
found extinction thresholds for dense cores found in Bolo- 
cam 1.1mm maps in the Perseus, Serpens, and Ophiuchus 
clouds of Ay ~8, 15, and 23, respectively (~120, 225, 
and 350 Mq pc 2 ), wit h a low probability of finding cores 
below these thresholds. lOnishi et alj (|1998f ) surveyed Tau- 
rus in C 18 and found a star forming column density 
threshold of 8xl0 21 cm~ 2 , which corresponds to a gas 
surfa ce density of 12 8.1 M pc 2 (Ay =8.5). Similarly, 
both I Johnstone et all (|2004[ ) and lAndre et al.1 (|2010l ) find 
t hresholds of A v ~10 (15 M^ p c 2 ) . 

iMouschovias k, Spitzerl (|1976l ) proposed the idea of a 
physical column density threshold corresponding to the 
central surface density above which the interstellar mag- 
netic field cannot support the gas fro m self gravita tional 
collapse. This was later modified by iMcKeg (|1989l) who 
consi dered the local ionization state s owing to UV radia- 
tion. IMouschovias fc Spitzerl (|1976f ) predicted that when 
clumps combine to form a large cloud complex, there ex- 
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ists a natural surface density threshold (£ C rit) for a given 
magnetic field: 



> 



80 



M Q pc 2 



D 



30 ^G 



(18) 



The total mean strength of the line-of-sight magnetic 
field (Bi os ) measure d in molecular c louds (n ~ 10 3 -10 4 



cm 3 ) is -10-20^G (|CrutcheH [19991 iTroland fc Crutchei 



2008). Since statistically Bi os « \B to t dHeiles fc Crutchei 



2005(), the corresponding total magnetic field, B to ti is ^20- 
40/xG. Using equation [IH the corresponding £ cr ;t >50- 
110 Mq pc 2 . A similar idea of a t hreshold at a par- 
ticular extinction was predicted by iMcKed (fl989) for 
photoionization-regulated star formation. This model pre- 
dicts that the rate of star formation is controlled by am- 
bipolar diffusion and therefore depends on the ionization 
levels in the cloud. Star formation in a "standard" ion- 
ization case will occur at Ay >4-8 mag, which trans- 
lates into a Scrit >60-120 M Q pc 2 . Both of these pre- 
dictions for a critical density of star formation are similar 
to E t h =129±14 Mq pc 2 . We note, however, that both 
these models are for parts of clouds that are in a "quasi- 
static" or turbulence-supported state. Alternatively, these 
parts of clouds may never become bound and are tran- 
sient ( Vazquez-Semadeni et al.l[2009T ). In this picture, the 
threshold would correspond to parts of molecular clouds 
that become gravitationally bound and form stars. 

3.3. Does the Lack of Resolution in Extragalactic Studies 
Explain the Discrepancy in Ssfr ? 

The third possibility is based on the fact that the ex- 
tragalactic relations are averaging over large scales which 
do not resolve the regions where stars are forming. Cur- 
rent "spatially resolved" extragala ctic measurements are 
still limited to scales of ~0.2-2 kpc dKennicutt et al.ll200"7t 
iBigiel et al.H2008HBlanc et al.ll2009ft : therefore, we cannot 
directly measure extragalactic SFRs on scales of galactic 
star forming regions. In any given spatially resolved ex- 
tragalactic measurement of Ssfr and S gas , the beam will 
contain a fraction of diffuse gas that does not trace star for- 
mation and a fraction of dense, star forming gas (/dense)- 
As discussed in Section |3~T1 the dense gas that is forming 
stars is not well traced by CO so there will be an excess of 
non-star forming gas in each beam measurement. A local 
example of this i s a stu dy of the Taurus molecular cloud; 
IGoldsmith etal] (12008ft found that 33% of diffuse 12 CO 
is contained in regions not associated with 13 CO, which 
is a more reliable tracer of dense, star forming gas (Sec- 
tion 13. ip . Even in the regions with Ay > 2 studied by the 
c2d and GB projects, star formation is highly concentrated 
to regions of high extinction (e.g., Figure[5]). Extragalactic 
studies averaging over hundreds of pc scales would include 
this diffuse gas, causing an increase in the amount of CO 
flux that is being counted as star forming gas. In order to 
better understand approximately how much gas is form- 
ing stars at present, a measurement of the fraction of gas 
that contains YSOs over a larger area on kpc scales in the 
G alaxy is neede d. 

lLada fe Ladal {2003) proposed the idea that clusters of 
stars which form in clumps located in giant molecular 
clouds are the fundamental building blocks of galaxies. 
The rate at which these stellar clusters form is set by the 



mechanisms that enable these clumps to condense out of 
th eir low den s ity pa rent cloud. A similar idea was explored 
by IWu et al.l (|2005f ). who proposed that there is a basic 
unit of clustered star formation with the following typical 
properties: Lin > 10 5 £©, i?donsc ~0.5 pc, and M vu . ~300- 
1000 M Q . As more of these basic units are contained in 
a galaxy, the SFR increases linearly. This linear correla- 
tion between SFRir and the mass of dense gas (Md en se) 
from HCN J=l-0 was seen in a sample of both spiral 
and lumin ous or ultra-luminous IR gala xies (LIRGS and 
ULIRGS) (|Gao fc Solomon! I2004bh . and IWu et all ((20051) 
showed that the same relation fit the Galactic massive 
dense clumps. It is therefore the dense gas tracers, such 
as HCN, that directly probe the volume of gas from which 
stars form in dense clumps and produce the star formation 
in external galaxies. 

If a linear relation between dense gas and the SFR is 
assumed at all S gas and Ssfr densities, how can we ex- 
plain the non-linear behavior of the Kennicutt-Schmidt 
SFR-gas relation? Let us suppose that the underlying 
star formation relation is what we actually observe in re- 
gions forming massive stars: a threshold around 129 M Q 
pc -2 and a roughly linear relation between dense gas and 
star formation above that threshold: 



J SFR 



oc E d 



(19) 



Let us also assume, the Kennicutt-Schmidt relation for 
the gas surface density averaged over large scales ((E gas )): 

^SFR OC \^Jgas ) 1 ' 4 . (20) 

Then, the fraction of gas above the threshold (/dense) 
would have to scale with mean surface density of all gas: 



— ^densc/ (^gas) 

oc S S fr/(S gas / 



(21) 



oc (So 



, 0.4 



Near the S t h of ~129 Mq pc~ 2 , the average Ssfr mea- 
sured on small scales is about ^3 Mq yr _1 kpc -2 (taking 
the average between Class I and Flat SED sources) and 

/dense IS about 1/40. When (Sgas) ~ 300£ t h, /dense ~ 1. 

At this point, all the gas is dense enough to form stars 
and star formation is most efficient, creating a maximal 
starb urst. Thi s is also where t he dense gas (|Wu et al.l 
2005|) and CO (|Kennicuttl Il998ah relations cross (see Fig- 
ure [5]). Above ^300£ t h, the only way to increase star 
formation efficiency is to make it more efficient even in 
the dense gas. This is possible because even in dense gas, 
the star formation rate per free-fall time is less than unity 
(jZuckerman fc Evansl[T97l iKrumholz fc Tanll2007f) . 

Is there local evidence for a preponderance of gas be- 
low the threshold? Complete maps of CO for the local 
0.5 kpc are not readily available, but we can measure the 
mass below and the mass above the threshold of Ay = 8 
for the 16 clouds with Spitzer coverage down to Ay = 2. 
The ratio of total mass lying below the threshold to the 
total mass above the threshold is 4.6. The massive star 
forming region, Orion, also has a similar ratio of 5.1 for 
the mass below over the mass above Ay = 8 (M. Heyer, 
unpublished data). A few clouds have been mapped to 
still lower levels and a factor of two more ma ss is found 
in Taurus, for example IGoldsmith et al.l (|2008f) . Alterna- 
tively, if we assume Orion to contain the largest reservoir 



Galactic Star Formation Rate and Gas Surface Density Relations 



11 



of molecular material within 0.5 kpc, we can derive the ra- 
tio of mass from 12 CO and Ay maps to get 6.4 (M. Heyer, 
unpublished data). Taken together, these factors make 
it plausible that there is 10 times more molecular mass 
than mass above the threshold. Furthermore, most of the 
gas within 0.5 kpc is atomic. If that is included, the pre- 
dictions of the extragalactic Kennicutt-Schmidt relation 
for total gas agree reasonabl y with the lo cal star forma- 
tion rate surface density (see Evansl (|2008f) and references 
therein). 

Finally, if the underlying star formation law in the dense 
gas is linear, then the arguments invoking the density de- 
pendence of the free-fall time to get a power of 1.5 (see 
Section [1]) are specious. In fact, the idea of a single free- 
fall time for a molecular cloud with an enormous range of 
densities is highly dubious in the first place. 

4. SUMMARY 

We investigate the relation between star formation rate 
(SFR) and gas surface densities in a sample of young stel- 
lar objects (YSOs) and massive dense clumps. Our YSO 
sample consists of objects located in 20 large molecular 
clouds from the Spitzer cores to disks (c2d) and Gould's 
Belt (GB) surveys. We estimate the £ gas in the c2d and 
GB clouds from Ay maps and Ssfr from the number of 
YSOs, assuming a mean mass and star formation timescale 
for each source. We also divide the clouds into evenly 
spaced contour levels of Ay. In each contour interval, we 
measure the £ gas .con and estimate the EgpR by counting 
only Class I and Flat SED YSOs which have not yet mi- 
grated from their birthplace. We use 12 CO and 13 CO gas 
maps of the Perseus and Ophiuchus clouds from the COM- 
PLETE survey to estimate £h 2 densities and compare to 
measurements from Ay maps. We also compare the c2d 
and GB low-mass star formin g regions t o a sam ple of mas- 
sive star forming clumps from lWu et al.l (|2010f ) . We derive 
SFRs from the total IR luminosity and use HCN gas maps 
to estimate SFR surface densities (£sFR m ) and gas surface 
densities (Ehcn) for the massive clumps. Our results are 
as follows: 

• The c2d and GB clouds lie above the extragalac- 
tic SFR-gas relations (e.g., Kennicutt-Schmidt 
Law) by factors of 9-17. We compare the to- 
tal cloud points to t he theoretical prediction of 
iKrumholz et alj ((2009) for galactic metallicity and 
a clumping factor of 1, corresponding to scales of 
100 pc, and find the clouds to lie above this predic- 
tion by a factor of —40. 

• We perform a follow up survey of suspicious YSOs 
(MISFITS) at the CSO using the HCO+ J=3-2 line 
transition as a dense gas tracer. We choose the 
youngest YSOs (Class I and Flat SED) that have 
not yet had time to migrate from their birthplace. 
These sources are spatially positioned at low ex- 
tinction levels, most are not clustered, and most lie 
outside the Ay peaks. In this paper, we present 
results for a total of 98 sources, including 45 Flat 
SED and 53 Class I YSOs (detailed results from 
the full survey will be published in a later paper). 
We find that 74% or 73 out of the 98 MISFITS 
observed to date are not detected in HCO + which 



indicates that they do not have a dense envelope of 
gas, and could be either later class YSOs or back- 
ground galaxies. These are a small fraction of the 
total number of YSOs in the sample, but they could 
bias the statistics upward at low E gas . 

• We divide the c2d and GB clouds into contours us- 
ing evenly spaced intervals of Ay fSection l2.2[) . We 
count only the youngest YSOs, removing any Class 
I or Flat SED YSOs (MI SFITS ) that are not de- 
tected in HCO+ ( Sections |2~XT1 and l2~l) . We find 
that the observed e xtrag alactic relations (|Kennicuttl 
1998b: iBigiel et al.ll2008l) under-predict the average 
E SF r of -9.7 Mq yr" 1 kpc" 2 by factors of -21-54 
and that our data lie a bove the theoretical relation 
(|Krumholz et a l. 2009) by —2 orders of magnitude. 

• We compare £ gas calculated from Ay maps to £h 2 
estimated from 12 CO in Section |3~T1 We find that 
the mass estimated from 12 CO may underestimate 
the S gas at E gas >200 M pc" 2 by >30% (Fig- 
ure [SJ. If the Eh 2 from 12 CO underestimates the 
H2 mass at S gas >200 Mq pc~ 2 , then this would 
effectively shift the extragalactic observed data to 
the rig ht above this thres hold, flattening the slope 
in the iKennicuttl (|1998bf ) relation. However, this 
small change is not enough to account for the dis- 
crepancy between Galactic and extragalactic mea- 
surements. 

• We also compare £r- 2 from 13 CO maps to S gas from 
Ay maps. If 13 CO traces the mass we find us- 
ing extinction maps, we would expect the ratio of 
13 CO/ 'Ay mass to be of order unity. However, we 
find the mass estimated from 13 CO to be lower than 
the Ay mass by factors of —4-5 (Section |3~T|) . 

• We find a steep decrease in SsFR/S gas (Figure [8]) 
and denote this as a star formation threshold (Sth)- 
In order to determine Sth, we fit a single power law 
and broken power law models to data for Class I 
and Flat SED YSOs and massive clumps. We find 
the best fit to the SFR-gas relation between YSOs 
and clumps to be a broken power law (Section l3.2.1l) 
with a break £ t h = 129±14 M Q pc -2 . We find a 
steep relation at S gas < E t h (slope of —4.6) and a 
linear relation at S gas > S t h with a slope of —1.1 
(Section GEO}. 

• Since the c2d and GB clouds are forming low-mass 
stars, and extragalactic studies are only able to use 
tracers that measure the light coming from mas- 
sive stars, the two star forming regimes might be- 
have differently, accounting for the large difference 
we measure. However, we find that both high and 
low-mass star forming regions in the Galaxy follow 
roughly the same linear relation above £ t h (Sec- 
tion l3~2"j). 

• A contributing factor to the difference seen be- 
tween Milky Way clouds and extragalactic measure- 
ments both on disk-averaged and spatially resolved 
scales is that extragalactic measurements average 
over large scales. These measurements include both 
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star forming gas and gas that is not dense enough 
to form stars. 



Assuming the Kennicutt-Schmidt relation and that 
the fundamental correlation between Esfr and the 
dense gas (Sdense) is linear, then the fraction of 
dense star- forming gas is proportional to (S gas ) ' 4 . 
When (Sgas) reaches ^300S t h, the fraction of dense 
gas is ~1, creating a maximal starburst. 
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Table 1 

Measured Quantities for Clouds. 



Cloud NVsOs.tot JVysOs.i W y sOs,f Distance SI A cloud Afgas, cloud £ ga3 , cloud SFR Ssfr 

(pc) (dog 2 ) (pc 2 ) (M ) (M (M (M yr" 

pc~ 2 ) Myr -1 ) kpc -2 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



Cha II 


24 





2 


178±18 


1.0 


9.9±2.0 


637.4±296.3 


64.3±26 


6.0±3.2 


0.61±0.35 


Lup I 


13 


2 


1 


150±20 


1.3 


8.9±2.4 


512.5±308.2 


57.9±31 


3.2±1.8 


0.37±0.22 


Lup III 


68 


2 


6 


200±20 


1.3 


15.4±3.1 


912.1±517.3 


59.1±31 


17.0±9.2 


1.10±0.63 


Lup IV 


12 


1 





150±20 


0.4 


2.5±0.7 


189.3±95.5 


75.1±32 


3. Oil. 6 


1.19±0.72 


Oph 


290 


27 


44 


125±25 


6.2 


29.6±11.8 


3115.3±1754.3 


105.4±41 


72.5±39.0 


2.45±1.65 


Per 


385 


76 


35 


250±50 


3.8 


73.2±29.3 


6585.3±3557.1 


90.0±32 


96.2±51.8 


1.32±0.88 


Ser 


224 


31 


21 


260±10 


0.8 


17. Oil. 3 


2336.7±640.2 


137.3±36 


56.0±30.2 


3.29±1.79 


AurN 


2 


1 





300±30 


0.1 


2.4±0.5 


223.9±51.9 


92.8±10 


0.5±0.3 


0.21±0.12 


Aur 


171 


43 


24 


300±30 


1.8 


50.0±10.0 


4617.5±1072.7 


92.4±10 


42.8±23.0 


0.86±0.49 


Ccp 


118 


30 


10 


300±30 


1.4 


38.0±7.6 


2610.3±168.5 


68.7±17 


29.5±15.9 


0.78±0.45 


Cha III 


4 


1 





200±20 


2.3 


28.0±5.6 


1326.0±386.2 


47.4±10 


1.0±0.5 


0.04±0.02 


Cha I 


89 


10 


12 


200±20 


0.8 


9.4±1.9 


857.3±206.3 


91.1±12 


22.2±12.0 


2.36±1.36 


CrA 


41 


7 


3 


130±25 


0.6 


3. Oil. 2 


279.2±114.0 


92.3±12 


10.2±5.5 


3.39±2.24 


IC5146E 


93 


13 


9 


950±80 


0.2 


61.4±10.3 


3365.2±872.9 


54.8±10 


23.2±12.5 


0.38±0.21 


IC5146NW 


38 


16 


3 


950±80 


0.3 


87.6±14.8 


5178.1±1257.3 


59.1±10 


9.5±5.1 


0.11±0.06 


Lup VI 


45 





1 


150±20 


1.0 


6.7±1.8 


454.9±141.4 


67.5±10 


11.2±6.1 


1.67±1.00 


Lup V 


43 








150±20 


1.7 


11.7±3.1 


704.7±223.5 


60.5±10 


10.8±5.8 


0.92±0.55 


Mus 


12 


1 





160±20 


0.9 


6.8±1.7 


335.1±109.1 


49.1±10 


3. Oil. 6 


0.44±0.26 


Sco 


10 


2 


1 


130±15 


1.4 


7.3±1.7 


620.6±17.4 


85.2±22 


2.5±1.3 


0.34±0.20 


Scr-Aqu 


1440 


146 


96 


260±10 


8.7 


179.5±13.8 


24441.3±3025.2 


136.2±13 


360.0±193.9 


2.01±1.09 


Cloud Averages 


156.1±71.5 


20.5±7.9 


13.4±5.2 


274.6±53.3 


1.8±0.5 


32.4±9.6 


2965.1±1204.9 


79.3±5.8 


39.0±17.9 


1.2±0.2 


Cloud Total 


3122.0 


409.0 


268.0 




36.0 


648.3 


59302.7 


91.5 


780.5 


1.2 


Data from Literature: 






















Taurus 


148 






137 


44 


252 


27207 


108 


37 


0.147 



Note. — Columns arc : (1) Cloud name; (2) Total number of YSOs at all Ay.; (3) Number of Class I objects at all Av ■', (4) Number of Flat SED objects at all 
Ay.; (5) Distances to each cloud.; (6) Solid angle.; (7) Area (pc 2 ).; (8) Mass (M©).; (9) Surface gas density (M© pc -2 ).; (10) Star formation rate (M© Myr" 1 ).; (11) 
Star formation rate density (M© yr -1 kpc~ 2 ).; (I) Total Av mass from Pineda ct al. (2010!) and YSO data from Rcbull ct 1^ | |2U1U| ). 



Table 2 

Measured Quantities for Clouds in Av Contours 



4^ 



Cloud -/VysOb.i -/VysOs,f Contour SI A con M con S gas ,con SFR, I SFR, F Esfr.i Ssfr.f 

levels 1 (dog 2 ) (pc 2 ) (M ) (M § (M (M (M q yr" 1 (M yr" : 

(mag) pc 2 ) Myr x ) Myr x ) kpc 2 ) kpc 2 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



Cha II 








5.2 


0.8 


7.78±1.57 


416.9±223.1 


53.6±26.6 


0.9±1.0* 


1.4±1.0* 


O.lzbO.l* 


r\ it 1 ft i "M 

U.2±u. 1 




Cha II 





1 


8.2 


0.2 


1.7o±U.3b 


lb2.o±oo.7 


92.3±27.b 


ft CI 1 1 i "i :■: 

U.9±1.U 


1.4il.O 


n C 1 ft c sfc 

u.o±u.b 


0.8±0.b 




Cha II 


f) 


1 


11.8 


0.03 


0.30±0.06 


a o ft 1 1 O T 

43.9±12.7 


1 A T ft 1 *~>ft A 

147.0±30.4 


ft ft 1 i n^k 

0.9±1.0 ¥ 


1.4±1.0 


O ft 1 o o 

3.0±3.3 


4.7±3.3 




Cha II 


1 





16.0 


0.01 


0.07±0.01 


14.0±3.8 


19o.9±o4.7 


0.9±1.0 


1.4±1.0 


i n c _l_ too 

12.6±13.8 


i n o_L_ i o o * 

19.2±13.8 




Lup I 


1 





6.0 


1.2 


8.05±2.15 


417.5±273.2 


51.9±31.0 


0.9±1.0 


1.4±1.0* 


0.1±0.1 


0.2±0.1* 




Lup I 








10.0 


0.1 


0.71±0.19 


77.8±31.0 


109.3±32.4 


0.9±1.0* 


1.4±1.0* 


1.3±1.4* 


2. 0±1. 4* 




Lup I 


f) 





16.0 


0.01 


0.09±0.02 


17.2±5.8 


184.1±39.0 


0.9±1.0* 


1.4±1.0* 


9.8±10.7* 


1 A ft 1 1 ft Tsfc 

14.9±10.7 




Lup III 


1 


2 


8.0 


1.2 


14.89±2.98 


815.8±490.4 


54.8±31.1 


0.9±1.0 


2.8±1.5 


O.liO.l 


0.2±1.4 




Lup III 





1 


14.0 


0.03 


0.42±0.08 


65.0±20.3 


153.1±36.6 


0.9±1.0* 


1.4±1.0 


2.1±2.4* 


3.3±2.4 




Lup III 


1 





20.0 


0.01 


0.13±0.03 


31.3±8.7 


249.0±48.0 


0.9±1.0 


1.4±1.0* 


7.2±8.0 


11.1±8.0* 




Lup IV 








8.0 


0.3 


2.26±0.60 


139.0±79.8 


61.4±31.2 


0.9±1.0* 


1.4±1.0* 


0.4±0.4* 


0.6±0.4* 




Lup IV 








14.0 


0.02 


0.17±0.04 


26.0±9.2 


156.6±36.5 


0.9±1.0* 


1.4±1.0* 


5.5±6.0* 


8.4±6.0* 




Lup IV 








23.0 


0.01 


0.09±0.02 


24.3±7.8 


266.8±46.9 


0.9±1.0* 


1.4±1.0* 


10.0±11.0* 


15.3±11.0* 




Oph 





2 


10.5 


5.6 


26.61±10.64 


2323.2±1421.2 


87.3±40.4 


0.9±1.0* 


2.8±1.5 


0.03±0.04* 


0.1±1.4 




Oph 


1 


3 


18.0 


0.4 


1.86±0.75 


368.1±170.7 


197.4±46.3 


0.9±1.0 


4.2±2.2 


0.5±0.5 


2.2±1.7 




Oph 


5 


5 


25.5 


0.1 


0.57±0.23 


182.3±80.1 


319.6±58.1 


4.5±2.5 


6.9±3.7 


8.0±5.4 


12.2±8.2 




Oph 


9 


14 


33.0 


0.1 


0.34±0.14 


147.5±63.8 


436.7±72.0 


8.2±4.5 


19.4±10.5 


24.2±16.4 


57.6±38.6 




Oph 


10 


12 


41.0 


0.04 


0.17±0.07 


94.3±40.8 


541.8±89.8 


9.1±5.0 


16.7±9.0 


52.2±35.4 


95.8±64.2 




Per 


3 





6.5 


2.7 


52.28±20.91 


3504.2±2158.9 


67.0±31.4 


2.7±1.5 


1.4±1.0* 


0.1±1.7 


0.03±0.02* 




Per 


15 


4 


11.0 


0.8 


15.39±6.16 


1880.9±910.5 


122.2±33.3 


13.6±7.5 


5.6±3.0 


0.9±3.9 


0.4±2.0 




Per 


18 


14 


15.5 


0.2 


3.79±1.52 


734.4±331.1 


193.6±40.3 


16.4±9.0 


19.4±10.5 


4.3±4.2 


5.1±3.7 




Per 


28 


11 


20.0 


0.1 


1.49±0.60 


388.8±170.4 


260.3±46.7 


25.5±13.9 


15.3±8.2 


17.0±11.6 


10.2±6.9 




Per 


8 


1 


24.5 


0.01 


0.22±0.09 


69.8±30.2 


316.9±52.8 


7.3±4.0 


1.4±1.0 


33.0±22.4 


6.3±4.5 


p' 

o 


Per 








30.0 


0.001 


0.02±0.01 


7.2±3.2 


403.0±81.0 


0.9±1.0* 


1.4±1.0* 


50.9±56.0* 


77.8±56.0* 


Ser 


2 


1 


10.2 


0.6 


13.34±1.03 


1588.5±479.4 


119.1±34.7 


1.8±1.0 


1.4±1.0 


0.1±1.4 


O.liO.l 


a 
p 
p: 


Ser 


2 


2 


14.5 


0.1 


2.54±0.20 


456.7±104.0 


179.9±38.5 


1.8±1.0 


2.8±1.5 


0.7±1.4 


1.1±1.4 


Ser 


7 


5 


18.8 


0.04 


0.91±0.07 


221.5±44.3 


243.2±44.9 


6.4±3.5 


6.9±3.7 


7.0±3.9 


7.6±4.1 


P 

a> 


Ser 


20 


9 


23.0 


0.01 


0.23±0.02 


70.1±13.3 


306.6±53.3 


18.2±9.9 


12.5±6.7 


79.6±44.0 


54.7±29.8 


Aur 


19 


6 


8.8 


1.8 


48.05±9.61 


4293.9±1001.2 


89.4±10.7 


17.3±9.4 


8.3±4.5 


0.4±4.4 


0.2±2.4 




Aur 


20 


15 


15.5 


0.1 


1.77±0.35 


274.0±60.1 


154.9±13.9 


18.2±9.9 


20.8±11.2 


10.3±6.0 


11.8±6.8 




Aur 


2 


1 


22.2 


0.01 


tt i o 1 ft ft A 

0.18±0.04 


A ft "7 1 1 1 O 

49.7±11.8 


it o it o ! or o 

282.8±35.8 


1.8±1.0 


1.4±1.0 


10.4±6.0 


7.9±5.7 




AurN 








5.2 


0.01 


0.46±0.09 


30.9±7.8 


67.4±10.4 


0.9±1.0* 


1.4±1.0* 


2.0±2.2* 


3.0±2.2* 




AurN 


1 





8.3 


0.1 


1.95±0.39 


193.0±44.2 


98.8±11.0 


0.9±1.0 


1.4±1.0* 


0.5±0.5 


0.7±0.5* 




Cep 


13 


3 


7.5 


1.3 


35.77±35.77 


2327.7±597.5 


331.9±52.1 


11.8±7.5 


4.2±3.0 


0.3±0.2 


O.liO.l 




Cep 


13 


3 


13.0 


0.1 


2.21±2.21 


282.6±63.3 


r nT it 1 rr i o 

507.2±51.3 


11.8±7.0 


4.2±3.0 


5.4±3.3 


1.9±1.4 




Cha I 


3 


6 


8.0 


0.6 


T A A _l_ 1 An 

7.44±1.4y 


c A A i _l_ ioc n 


to i _l_ ii n 


2.7±1.5 


8.3±4.5 


. 4 ± 1 . 7 


l.li2.4 




Cha I 


7 


5 


14.0 


0.1 


i T rr _l_ a or 

1.75±u.3o 


it C r* r\ _|_ C T C 

25o.U±57.5 


~\ A C C _l_ 1 C ft 

14b.b±I5.U 


6.4±3.5 


6.9±3.7 


3.6±2.6 


4.0i2.3 




Cha I 





1 


21.0 


0.02 


0.22±0.04 


56.0±13.1 


253.6±30.9 


0.9±1.0* 


1.4±1.0 


4.1±4.5* 


6.3±4.5 




Cha III 


1 





5.0 


2.2 


26.85±5.37 


1228.5±364.1 


45.8±10.0 


0.9±1.0 


1.4±1.0* 


0.03±0.04 


0.05±0.04* 




Pha TTT 








8.0 


0.1 




Q7 n 


or 9-Li n a 


9il 0* 


1.4±1.0* 


0.8±0.9* 


1.2±0.9* 




CrA 


2 


1 


9.3 


0.5 


2.45±0.94 


162.8±68.1 


66.5±11.0 


1.8±1.0 


1.4±1.0 


0.7±1.4 


0.6±0.4 




CrA 


1 


2 


16.7 


0.09 


0.47±0.18 


87.2±34.7 


183.5±18.2 


0.9±1.0 


2.8±1.5 


1.9±2.1 


5.8±3.9 




CrA 


4 





24.0 


0.02 


0.10±0.04 


29.1±11.6 


289.8±28.4 


3.6±2.0 


1.4±1.0* 


36.2±24.2 


13.8±9.9* 




IC5146E 








4.7 


0.2 


50.60±8.52 


2424.3±673.2 


47.9±10.6 


0.9±1.0* 


1.4±1.0* 


0.02±0.02* 


0.03±0.02* 




IC5146E 


11 


6 


7.4 


0.04 


10.82±1.82 


940.8±204.1 


86.9±11.9 


10.0±5.5 


8.3±4.5 


0.9±3.3 


0.8±2.4 




IC5146NW 


7 





5.5 


0.3 


73.75±12.42 


3831.8±989.0 


52.0±10.2 


6.4±3.5 


1.4±1.0* 


0.1±2.6 


0.02±0.01* 




IC5146NW 


8 


3 


9.0 


0.05 


13.83±2.33 


1346.3±275.4 


97.3±11.3 


7.3±4.0 


4.2±2.2 


0.5±2.8 


0.3±1.7 




Lup V 








4.5 


1.3 


9.08±2.42 


514.5±166.2 


56.7±10.3 


0.9±1.0* 


1.4±1.0* 


O.liO.l* 


0.2±0.1* 




Lup V 








7.0 


0.4 


2.57±0.69 


190.0±57.5 


73.9±10.6 


0.9±1.0* 


1.4±1.0* 


0.4±0.4* 


0.5±0.4* 




Lup VI 








4.5 


0.6 


3.83±1.02 


237.4±75.4 


62.0±10.7 


0.9±1.0* 


1.4±1.0* 


0.2±0.3* 


0.4±0.3* 




Lup VI 








9.0 


0.4 


2.90±0.77 


216.4±65.8 


74.6±10.9 


0.9±1.0* 


1.4±1.0* 


0.3±0.3* 


0.5±0.3* 




Mus 


1 





4.5 


0.8 


5.86±1.47 


259.2±88.0 


44.2±10.2 


0.9±1.0 


1.4±1.0* 


0.2±0.2 


0.2±0.2* 




Mus 








7.0 


0.1 


0.95±0.24 


75.4±21.4 


79.1±10.7 


0.9±1.0* 


1.4±1.0* 


1.0±1.0* 


1.5±1.0* 




Sco 


1 


1 


7.5 


1.2 


6.35±6.35 


484.7±130.6 


456.4±63.9 


0.9±1.0 


1.4±1.0 


0.1±0.2 


0.2±0.2 




Sco 


1 





17.0 


0.2 


0.94±0.94 


135.8±34.3 


567.1±57.6 


0.9±1.0 


1.4±1.0* 


l.Oil.l 


1.5±1.1* 




Ser-Aqu 








7.0 


3.0 


61.01±4.69 


5509.9±782.3 


90.3±10.8 


0.9±1.0* 


1.4±1.0* 


0.02±0.02* 


0.02±0.02* 




Ser-Aqu 


9 


4 


12.0 


4.6 


95.23±7.33 


12818.2±1540.2 


134.6±12.4 


8.2±4.5 


5.6±3.0 


0.1±3.0 


0.1±2.0 





Table 2 — Continued 



Cloud 
(1) 


NySOs,I 

(2) 


NySOs,F 

(3) 


Contour 
levels 1 
(mag) 
(4) 


a 

(dog 2 ) 
(5) 


A con 
(pc 2 ) 

(6) 


M con 
(Mq) 

(7) 


s gas , con 
(Mo 
pc" 2 ) 
(8) 


SFR, I 

(M e 

Myr" 1 ) 
(9) 


SFR, F 
(Me 

Myr- 1 ) 
(10) 


(M yr" 1 
kpc" 2 ) 
(11) 


^SFR.F 

(M e yr" 1 
kpc" 2 ) 
(12) 


Ser-Aqu 


16 


12 


17.0 


0.6 


12.57±0.97 


2658.5±298.5 


211.4±17.3 


14.5±8.0 


16.7±9.0 


1.2±4.0 


1.3±3.5 


Scr-Aqu 


31 


20 


22.0 


0.3 


6.66±0.51 


1929.6±219.8 


289.6±24.3 


28.2±15.4 


27.8±15.0 


4.2±5.6 


4.2±4.5 


Scr-Aqu 


31 


25 


27.0 


0.1 


2.88±0.22 


1036.8±125.3 


359.6±33.5 


28.2±15.4 


34.7±18.7 


9.8±5.6 


12.0±6.6 


Scr-Aqu 


50 


24 


33.0 


0.1 


1.10±0.08 


488.4±65.1 


442.5±48.2 


45.5±24.9 


33.3±18.0 


41.2±22.8 


30.2±16.4 


Contour Averages 


6.6±1.2 


4.1±0.7 


-999.0±-999.0 


0.6±0.1 


10.6±2.5 


972.1±248.0 


188.1±17.7 


6.0±1.1 


5.7±1.0 


7.7±2.0 


8.6±2.3 



Note. — Columns arc : (1) Cloud name.; (2) Number of Class I YSOs in contour level.; (3) Number of Flat SED YSOs in contour level.; (4) Ay contour level in mag at 
which mass measurement was made. The contour levels start at As/— 2 or the cloud completeness limit and increase in even intervals to the listed contour level.; (5) Solid 
angle.; (6) Area in contour level (pc -2 ).; (7) Mass in contour level (Mq),; (8) Surface gas density in contour level (Mq pc -2 ).; (9) Star formation rate in contour level (Mq 
Myr -1 ) for Class I YSOs. Asterisks denote that measurement is an upper limit.; (10) Star formation rate in contour level (Mq Myr -1 ) for Flat SED YSOs. Asterisks denote 
that measurement is an upper limit.; (11) Star formation rate density in contour level (Mq yr - 1 kpc -2 ) for Class I YSOs. Asterisks denote that measurement is an upper 
limit.; (12) Star formation rate density in contour level (Mq (yr -1 kpc -2 ) for Flat SED YSOs. Asterisks denote that measurement is an upper limit.; (I) Cont our levels start 
at Ay — 2 for all clouds except for Serpens and Ophiuchus which are covered by the c2d survey completely down to Av — 6 and 3 as discussed in Scction |2.1| 
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Table 3 

Properties of Suspicious YSOs and MISFITS 



Cloud 


RA 


DEC 


Classification 


a 


/ T MB dV 




SED 


A v 


comments 




J2000 


.T2000 






(K km s" 1 ) 


(K) 


class 


(mag) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0) 


(10) 



Aur 


04 


18 


21 


27 


+38 


01 


35 


88 


YSOc 






-0 


09 




<0 


23 


Flat 


4 


Aur 


04 


19 


44 


67 


+38 


11 


21 


98 


YSOc.star+dust(IRl) 


-0 


07 




<0 


23 


Flat 


7 


Aur 


04 


29 


40 


02 


+35 


21 


08 


95 


YSOc 


_star+dust(IRl) 





51 




<0 


31 


I 


8 


Aur 


04 


30 


14 


96 


+36 


00 


08 


53 


YSOc 


.red 




1 


77 


0.42+0.09 


0.69+0 


15 


I 


8 


Aur 


04 


30 


23 


83 


+35 


21 


12 


35 


YSOc 


.red 







61 




<0 


30 


I 


8 


Aur 


04 


30 


41 


17 


+35 


29 


41 


08 


YSOc 


.red 




1 


49 


0.76+0.07 


1.49+0 


14 


I 


7 


Aur 


04 


30 


44 


23 


+35 


59 


51 


16 


YSOc 






1 


08 


0.76+0.15 


0.78+0 


15 


I 


8 


Aur 


04 


30 


48 


52 


+35 


37 


53 


76 


YSOc 


.red 




1 


46 


1.01+0.13 


1.15+0 


15 


I 


7 


Aur 


04 


30 


56 


62 


+35 


30 


04 


55 


YSOc 


.red 




2 


35 


0.49+0.07 


1.00+0 


15 


I 


7 


Cop 


22 


29 


33 


35 


+75 


13 


10 


01 


YSOc 


.red 







20 




<0 


42 


Flat 


6 


Cep 


22 


35 


00 


82 


+75 


15 


30 


42 


YSOc 


_star+dust(IR2) 


-0 


29 




<0 


46 


Flat 


6 


Cep 


22 


35 


14 


09 


+75 


15 


02 


01 


YSOc 


.red 







36 




<0 


42 


I 


6 


Cep 


21 


01 


36 


07 


+68 


08 


22 


54 


YSOc 






-0 


21 




<0 


42 


I 


5 


Cep 


21 


01 


43 


89 


+68 


14 


03 


31 


YSOc 


.red 







14 




<0 


40 


I 


6 


Cep 


21 


02 


14 


00 


+68 


07 


30 


80 


YSOc 


.red 







49 




<0 


40 


I 


6 


Cep 


21 


02 


21 


22 


+67 


54 


20 


28 


YSOc 


.red 







68 


0.89+ 0.18 


1.74+ 


34 


I 


9 


Cep 


21 


02 


21 


22 


+67 


54 


20 


28 


YSOc 


.red 







68 


0.72+ 0.11 


2.33+ 


34 


I 


9 


Cep 


21 


02 


21 


30 


+68 


04 


30 


11 


YSOc 


.PAH 


— cm 





52 




<0 


42 


I 


5 


Cep 


21 


02 


59 


40 


+68 


06 


32 


24 


YSOc 


.red 







65 




<0 


42 


I 


5 


IC5146E 


21 


52 


46 


58 


+47 


12 


49 


32 


YSOc 


_star+dust(IR2) 


-0 


19 




<0 


34 


Flat 


5 


IC5146E 


21 


53 


36 


24 


+47 


10 


27 


84 


YSOc 


.star 4 


dust(IRl) 


-0 


12 




<0 


30 


Flat 


6 


IC5146E 


21 


54 


18 


70 


+47 


12 


09 


73 


YSOc 


_star+dust(IR2) 


-0 


23 




<0 


26 


Flat 


4 


IC5146E 


21 


52 


14 


36 


+47 


14 


54 


60 


YSOc.star+dust(IR2) 





67 




<0 


28 


I 


4 


IC5146E 


21 


52 


37 


78 


+47 


14 


38 


40 


YSOc_star+dust(IRl) 





64 


1.59+ 0.25 


1.80+ 


28 


I 


5 


IC5146E 


21 


53 


06 


94 


+47 


14 


34 


80 


YSOc 









34 


0.71+ 0.27 


0.66+ 


25 


I 


5 


IC5146E 


21 


53 


55 


70 


+47 


20 


30 


13 


YSOc 


.PAH 


—em 


1 


59 




<0 


34 


I 


4 


IC5146NW 


21 


45 


31 


22 


+47 


36 


21 


24 


YSOc 









13 


0.44+ 0.19 


0.61+ 


26 


Flat 


5 


IC5146NW 


21 


44 


43 


08 


+47 


46 


43 


08 


YSOc 


.red 







63 


1.92+ 0.11 


4.18+ 


23 


I 


4 


IC5146NW 


21 


44 


48 


31 


+47 


44 


59 


64 


YSOc 


.red 




1 


83 


0.65+ 0.05 


3.11+ 


23 


I 


5 


IC5146NW 


21 


44 


53 


98 


+47 


45 


43 


56 


YSOc 


_star+dust(IRl) 





76 


0.70+ 0.11 


1.75+ 


28 


I 


4 


IC5146NW 


21 


45 


02 


64 


+47 


33 


07 


56 


YSOc 


.red 




1 


18 


0.77+ 0.22 


1.03+ 


30 


I 


4 


IC5146NW 


21 


45 


08 


31 


+47 


33 


05 


77 


YSOc 


.red 







74 


3.72+ 0.53 


2.07+ 


30 


I 


4 


IC5146NW 


21 


45 


27 


80 


+47 


45 


50 


40 


YSOc 


.starH 


dust(IR4) 





42 




<0 


36 


I 


3 


IC5146NW 


21 


47 


06 


02 


+47 


39 


39 


24 


YSOc 


.red 





43 


0.40+ 0.14 


0.70+ 


25 


I 


5 


Lup 1 


15 


38 


48 


35 


-34 


40 


38 


24 


YSOc 


.PAH 


—cm 





31 




<0 


42 


I 


3 


Lup 1 


15 


43 


02 


29 


-34 


44 


06 


22 


YSOc 


.star 4 


dust(IRl) 





14 




<0 


38 


Flat 


<2 


Lup III 


16 


07 


03 


85 


-39 


11 


11 


59 


YSOc 


.star 4 


dust(IRl) 


-0 


14 




<0 


30 


Flat 


<2 


Lup III 


16 


07 


08 


57 


-39 


14 


07 


75 


YSOc 




-0 


01 




<0 


32 


Flat 


<2 


Lup III 


16 


07 


54 


73 


-39 


15 


44 


49 


YSOcrcd 




-0 


15 




<0 


28 


Flat 


2 


Lup IV 


16 


02 


21 


61 


-41 


40 


53 


70 


YSOcPAH 


em 





56 




<0 


36 


I 


4 


Lup VI 


16 


24 


51 


78 


-39 


56 


32 


66 


YSOc 









22 




<0 


48 


Flat 


8 


Oph 


16 


21 


38 


72 


-22 


53 


28 


26 


YSOc 


.star+dust(IRl) 





99 




<0 


35 


I 


<3 


Oph 


16 


23 


40 


00 


-23 


33 


37 


36 


YSOc 









01 




<0 


37 


Flat 


3 


Oph 


16 


44 


24 


27 


-24 


01 


24 


56 


YSOc 


.PAH 


— cm 





27 




<0 


35 


Flat 


<3 


Oph 


16 


45 


26 


65 


-24 


03 


05 


41 


YSOcred 







37 




<0 


36 


I 


<3 


Oph 


16 


21 


45 


13 


-23 


42 


31 


03 


YSOc.star+dust(IRl) 





30 




<0 


36 


I 


9 


Oph 


16 


31 


31 


24 


-24 


26 


27 


87 


YSOc 


_star+dust(IR4) 


-0 


24 




<0 


38 


Flat 


<3 


Oph 


16 


25 


27 


50 


-24 


36 


47 


55 


YSOc 


.star -f 


dust(IRl) 





06 




<0 


42 


Flat 


7 


Oph 


16 


23 


32 


22 


-24 


25 


53 


82 


YSOc 


_star+dust(IR2) 


-0 


04 


0.36+ 0.17 


0.66+ 


31 


Flat 


4 


Oph 


16 


22 


20 


99 


-23 


04 


02 


35 


YSOc 


.PAH 


— cm 





17 




<0 


44 


Flat 


4 


Oph 


16 


23 


05 


43 


-23 


02 


50 


73 


YSOc 


_star+dust(IR2) 


-0 


27 




<0 


46 


Flat 


4 


Oph 


16 


23 


06 


86 


-22 


57 


30 


01 


YSOc 






-0 


19 




<0 


44 


Flat 


5 


Oph 


16 


23 


40 


00 


-23 


33 


37 


36 


YSOc 









01 




<0 


46 


Flat 


4 


Per 


03 


25 


19 


52 


+30 


34 


24 


16 


YSOc 






-0 


11 




<0 


27 


Flat 


2 


Per 


03 


26 


37 


47 


+30 


15 


28 


08 


YSOc 


.red 







99 


0.20+0.04 


0.77+0 


14 


I 


2 


Per 


03 


26 


37 


47 


+30 


15 


28 


08 


YSOc 


.red 







99 


0.11+0.03 


0.49+0 


14 


I 


2 


Per 


03 


28 


34 


49 


+31 


00 


51 


10 


YSOc.star+dust(IRl) 





89 


0.31+0.05 


0.95±0 


15 


I 


6 


Per 


03 


28 


34 


94 


+30 


54 


54 


55 


YSOc 









01 




<0 


32 


Flat 


3 


Per 


03 


29 


06 


05 


+30 


30 


39 


19 


YSOc 


.red 







72 




<0 


30 


I 


2 


Per 


03 


29 


51 


82 


+31 


39 


00 


03 


red 






3 


34 


3.22±0.20 


2.41+0 


15 


I 


6 


Per 


03 


30 


22 


45 


+31 


32 


40 


53 


YSOc 


_star+dust(IR2) 





35 




<0 


34 


I 


3 


Per 


03 


30 


38 


21 


+30 


32 


11 


93 


YSOc.star+dust(IR2) 


-0 


10 




<0 


29 


Flat 


4 


Per 


03 


31 


14 


70 


+30 


49 


55 


40 


YSOc 


_star+dust(IRl) 


-0 


09 




<0 


30 


Flat 


2 


Per 


03 


31 


20 


98 


+30 


45 


30 


06 


YSOc 


.red 




1 


10 


2.70±0.17 


2.68+0 


17 


I 


5 


Per 


03 


44 


24 


84 


+32 


13 


48 


36 


YSOc 


.red 




1 


69 




<0 


31 


I 


6 


Per 


03 


44 


35 


34 


+32 


28 


37 


18 


YSOc 


.red 




-0 


09 




<0 


32 


Flat 


3 


Per 


03 


45 


13 


82 


+32 


12 


10 


00 


YSOcrcd 







43 




<0 


28 


I 


3 


Per 


03 


47 


05 


43 


+32 


43 


08 


53 


YSOc 


.red 







48 


0.93+0.09 


1.48+0 


14 


I 


5 


Sco 


16 


40 


58 


27 


-09 


35 


19 


76 


YSOc 


.red 







66 




<0 


44 


I 


12 


Sco 


16 


48 


28 


85 


-14 


14 


30 


45 


YSOc 


.PAH 


— cm 





48 




<0 


48 


I 


5 


Sco 


16 


22 


04 


35 


-19 


43 


20 


76 


YSOc 









02 




<0 


66 


Flat 


6 


Ser 


18 


28 


41 


87 


-00 


03 


21 


34 


YSOc.star+dust(IRl) 





14 




<0 


47 


Flat 


8 


Ser 


18 


28 


44 


78 


+00 


51 


25 


79 


YSOc 


.red 




1 


05 


0.86+0.23 


0.83±0 


22 


I 


8 


Ser 


18 


28 


44 


96 


+00 


52 


03 


54 


YSOc 


.red 




1 


27 


1.61+0.34 


1.36+0 


28 


I 


8 


Ser 


18 


29 


16 


18 


+00 


18 


22 


71 


YSOc 






-0 


13 


1.28+0.29 


0.98+0 


22 


Flat 


7 


Ser 


18 


29 


40 


20 


+00 


15 


13 


11 


YSOc 


.starH 


dust(IRl) 





68 




<0 


50 


I 


<6 


Ser 


18 


30 


05 


20 


+00 


41 


04 


58 


red 




1 


24 




<0 


44 


I 


<6 


Ser 


18 


28 


44 


01 


+00 


53 


37 


93 


YSOc 


.red 







29 




<0 


42 


Flat 


7 



sclf-rcvcrscd 



sclf-rcvcrscd 



double peak 
double peak 



double peak 
double peak 
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Table 3 — Continued 



Cloud 


RA 


DEC 


Classification 


a 


/ T MB dV 


2mb 


SED 


A v 


comments 




J2000 


J2000 






(K km s" 1 ) 


(K) 


class 


(mag) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



Ser 




18 


29:27.35 


+00 


38:49 


75 


YSOc 


0.24 


<0 


38 


Flat 


13 


Ser 




18 


29:31.96 


+01 


18:42 


91 


YSOc.star+dust(IRl) 


0.32 


0.63+ 0.27 0.69+ 


30 




11 


Ser- 


-Aqu 


18 


13:45.05 


-03 


26:02 


67 


YSOc_star+dust(IRl) 


0.39 


<0 


58 




6 


Ser- 


-Aqu 


18 


27:03.33 


-02 


45:33 


42 


YSOcrcd 


0.44 


<0 


54 




5 


Scr- 


-Aqu 


18 


29:16.80 


-01 


17:30 


68 


YSOc 


0.78 


<0 


54 




10 


Ser- 


-Aqu 


18 


30:32.48 


-03 


50:01 


21 


YSOc.star+dust(MPl) 


0.37 


<0 


54 




9 


Ser- 


-Aqu 


18 


33:03.49 


-02 


08:42 


53 


YSOc.PAH-cm 


1.25 


<0 


56 




8 


Scr- 


-Aqu 


18 


37:39.24 


-00 


25:35 


18 


YSOc.star+dust(IRl) 


0.72 


<0 


56 




7 


Ser- 


-Aqu 


18 


37:46.90 


-00 


01:55 


83 


YSOc.PAH-cm 


1.19 


<0 


52 




8 


Ser- 


-Aqu 


18 


37:52.77 


-00 


23:03 


10 


YSOc.star+dust(MPl) 


0.59 


<0 


50 




7 


Scr- 


-Aqu 


18 


37:55.79 


-00 


23:31 


59 


YSOc 


1.14 


<0 


58 




8 


Ser- 


-Aqu 


18 


05:31.11 


-04 


38:09 


63 


YSOc 


0.23 


<0 


60 


Flat 


9 


Ser- 


-Aqu 


18 


10:28.90 


-02 


37:42 


79 


YSOc 


-0.18 


<0 


62 


Flat 


6 


Scr- 


-Aqu 


18 


26:32.81 


-03 


46:27 


26 


YSOcred 


0.08 


<0 


60 


Flat 


10 


Ser- 


-Aqu 


18 


27:24.87 


-03 


58:21 


15 


YSOc 


-0.22 


<0 


60 


Flat 


10 


Scr- 


-Aqu 


18 


28:09.49 


-02 


26:31 


95 


YSOc.star+dust(IR2) 


-0.11 


<0 


58 


Flat 


5 


Scr- 


-Aqu 


18 


29:16.72 


-01 


17:36 


02 


YSOc_star+dust(MPl) 


-0.08 


<0 


62 


Flat 


10 


Ser- 


-Aqu 


18 


30:06.06 


-01 


10:19 


33 


YSOc 


-0.15 


<0 


56 


Flat 


6 


Scr- 


-Aqu 


18 


30:13.01 


-01 


25:36 


64 


YSOc.star+dust(IR2) 


-0.23 


<0 


60 


Flat 


8 


Ser- 


-Aqu 


18 


36:02.64 


-00 


02:20 


70 


YSOc_star+dust(MPl) 


-0.28 


<0 


58 


Flat 


8 


Scr- 


-Aqu 


18 


38:55.77 


-00 


23:40 


81 


YSOc 


-0.05 


<0 


64 


Flat 


7 


Scr- 


-Aqu 


18 


40:12.06 


+00 


29:27 


74 


YSOcrcd 


-0.07 


<0 


50 


Flat 


8 



Note. — Columns ar e : (1) Cloud (2) S ource Right Ascension in J2000 coordinates; (3) Source Declination in J2000 coordinates; (4) 
Source classification (sec Evans ct al. (2009)); (5) Spectral Index, extinction corrected values for c2d clouds only; (6) Integrated main beam 
HCO~*~line intensity.; (7) Main beam HCO~^line temperature, upper limits arc computed as 2a rms ] (8) SED class based on IGreene et ail 
(1994); (9) Av in source position. Values that are found outside the Ay map completeness limit are given as < limit; (10) Line profile 
comments. 



Table 4 

A v , 12 CO, and 13 CO Masses and E gas for Per and Oph Clouds 



Cloud 


M cIoud , gas 


M cloud.l2CO 


M cloud.l3CO 


^cloud,gas 


E cloud,12cO 


S cloud.l3C0 




(M ) 


(M ) 


(M s ) 


(Me 


(M e 


(M 










PC" 2 ) 


PC" 2 ) 


PC" 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Per 


5997+3387 


9657+2416 


1073+110 


82+33 


132+33 


15+2 


Oph 


2270+1533 


2596+659 


348+39 


77+42 


88+22 


12 + 1 


Data from Literature: 














Taurus 1 


27207 


16052 




108 


64 





Note. — (1) Cloud name.; (2) Mass from Av map (Mq) where there is positive 12 CO and 13 CO emission.; (3) 
12 CO Mass (M ).; (4) 13 CO Mass (Mq).; (5) Surface gas density from A v map (Mq pc~ 2 ).; (6) 12 CO Surface 
gas density (Me pc~ 2 ).; (7) 13 CO Surface gas density (Me pc~ 2 ).; (I) Combined 12 CO and 13 CO mass from 
IColdsmith et al.l 120081) and A v mass from HPineda et alJ|201Cf l. 
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Table 5 

Av, 12 CO, and 13 CO Masses and E gas for Per and Oph Clouds in Ay Contours 

Cloud Contour M co „, gaB M con i 2co M con i 3co Sco„,ga 3 £ con i 2co £ con ,i3 Cl 

levels (M ) (M ) (M ) (M (M (M 

(mag) pc 2 ) pc 2 ) pc 2 ) 

(1) (2) (3) (4) (5) (6) (7) (8) 



Per 


6.5 


3144.7±2068.5 


6050.8±1514.5 


518.0±54.0 


60.2±31 


115.7±28 


9.9±1 




11.0 


1774.0±875.5 


2555.9±639.3 


349.1±35.3 


115.3±33 


166.1±41 


22.7±2 




15.5 


652.9±302.6 


654.0±163.6 


120.3±12.1 


172.1±40 


172.4±43 


31.7±3 




20.0 


369.7±163.5 


343.7±86.0 


76.2±7.6 


247.6±46 


230.2±57 


51.0±5 




24.5 


49.5±23.0 


49.4±12.4 


8.9±0.9 


224.9±52 


224.5±56 


40.5±4 




30.0 


6.1±2.8 


3.3±0.8 


0.7±0.1 


339.3±80 


187.7±46 


37.6±3 


Oph 


10.5 


1577.2±1246.7 


1982.0±505 


183.3±23 


76.8±40 


74.5±19 


6.9±0 


18.0 


349.4±164.3 


328.4±82 


71.1±7 


187.4±46 


176.1±44 


38.1±3 




25.5 


161.3±72.5 


135.9±34 


42.6±4 


282.8±58 


238.3±59 


74.7±7 




33.0 


111.5±50.8 


97.1±24 


32.2±3 


330.3±72 


287.5±71 


95.3±9 




41.0 


70.6±32.3 


52.9±13 


18.6±1 


405. 9 ±89 


304.0±76 


106.6±10 





Note. — (1) Cloud name; (2) Av contour level in mag at which mass measurement was made. The contour 
levels start at Ay— 2 or the cloud completeness limit and increase in even intervals to the listed contour level.; (3) 
A v Mass (M ) where there is positive 12 CO and 13 CO emission.; (4) 12 CO Mass (M Q ).; (5) 13 CO Mass (M ).; 
(6) Av Surface gas density (Mq pc -2 ).; (7) 12 CO Surface gas density (Mq pc -2 ).; (8) 13 CO Surface gas density 
(M pc- 2 ). 



Table 6 
Massive Clumps HCNJ=(l-0) 

Source log Ehcn log Ssfr ir 

(Mq pc" 2 ) (M yr- 1 kpc" 2 ) 
(1) (2) (3) 



W3(OH) 


3.39+0, 


13 


1 


.35+0.12 


RCW142 


3.40+0, 


14 


1 


.08+0.13 


W28A2(1) 


3.66+0, 


14 


2 


.12+0.13 


G9. 62+0. 10 


3.28+0, 


13 


1 


,45+0.14 


G10. 60-0.40 


3.32+0. 


12 


1 


,83+0.12 


G12.21-0.10 


2.63+0, 


24 





.28+0.23 


G13.87+0.28 


2.28+0, 


15 





.61+0.13 


G23.95+0.16 


2.28+0, 


25 





,79+0.21 


W43S 


2.63+0, 


14 


1 


.51+0.14 


W44 


2.79+0, 


18 


1 


.00+0.16 


G35.58-0.03 


3.41+0, 


24 





,89+0.22 


G48. 61+0.02 


1.98+0, 


14 





.75+0.13 


W51M 


3.14+0. 


17 


1 


.53+0.17 


S87 


2.76+0, 


16 





.97+0.12 


S88B 


2.68+0, 


19 


1 


.31+0.15 


K3-50 


3.26+0, 


12 


1 


.75+0.13 


ONI 


2.83+0, 


14 





,73+0.13 


ON2S 


2.76+0, 


16 


1 


.44+0.14 


W75N 


3.13+0. 


13 


1 


.50+0.12 


DR21S 


3.20+0, 


13 


1 


,75+0.12 


W75(OH) 


3.21+0. 


13 





.62+0.13 


CEPA 


3.44+0. 


17 


2 


,00+0.13 


IRAS20126 


2.98+0, 


18 


1 


.22+0.13 


IRAS20220 


2.63+0, 


21 





.29+0.18 


IRAS23385 


2.79+0, 


19 





,43+0.15 


Clump Average 


3.12+0.16 


1.44+0.15 



Note. — Columns arc : (1) Source name (2) Surface 
gas density; (3) SFR surface gas density. 
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Fig. 1. — Ssfr is shown versus the S gas for c2d and GB clouds (cyan squares). All cloud S gas are measured above Ay > 2 (or the cloud 
com pleteness limit , see S ection 12.111 . We also include an estimate for the Tau rus molecular cloud (black square) which includes YSO counts 
fromlRebuII et al. l d2010 | ) and an Ay > 2 gas mass from IPineda et alj l|2010l) . Extragala ctic observed rel a tions are shown for the sample of 
Kennicutt ( 1998bf) and IBieiel et aLI 1 120081) (blue solid and red lines, respectively). The Krumholz et al. (2009) prediction for the total (H 
I+CO) gas star formation law for the galactic metallicity and a clumping factor of 1 corresponding to ~100 pc scales is also shown (orange 
line). 



20 



Heiderman et al. 



2 7 13 ^ ( mo< 3) 18 24 29 




: ° 

I I I I I I I I I I I I I L 
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RA (J2000) 

Fig. 2. — An example of the S gas measurement method in the Perseus molecular cloud from the c2d survey. The grayscale image is the 
extinction map with black contours ranging from 2-29 in intervals of 4.5 mag The yellow filled circles are Flat SED sources and the red filled 
circles are Class I sources. Sources that have an open star correspond to suspicious YSOs (MISFITS) that were observed in HCO+ J=3-2 at 
the CSO and were not detected. We measure the S gas from each map in each contour of extinction. Contours are spaced in intervals wider 
than the extinction map beam size of 270". To estimate SFR, we count the YSOs in corresponding contour levels (Section I2,2|l , 
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1.0 1.5 2.0 2.5 3.0 3.5 4.0 

"og 2 gos (M pc" 2 ) 

Fig. 3. — Gas surface densities measured from extinction maps and SFRs estimated from Class I (green stars) and Flat SED (magenta 
stars) YSO number counts in c2d and Gould's belt clouds are shown. For contour levels that do not contain any YSOs, we calculate an upper 
limi t for that region us ing one YSO (open inverted triangles). Extragalactic observed relations are shown for the sample of Kcnnicutt ( 1998b) 
and lBigiel et all (2008) (blue solid and red lines, respectively). 
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Fig. 4 . — top panels : 12 CO integrated intensity versus visual extinction (Ay ) for Per (left) and Oph (right). The standard Xco~f ac tor 
fit from [Bloemen et al.l i 19861) is shown by the dashed grey lines (Section 13. It . bottom panels: 13 CO column densities versus visual 
extinction (Ay) for Per (left) and Oph (right). The average H2-to- 13 CO abundance ratio from the literature is shown by the grey dashed 
lines (Section [XT} . 
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Fig. 5. — The ratio of H2 gas surface densities from CO compare to that estimated from Ay maps (S gaa ). The cyan squares and circles are 
points for the Oph and Per clouds, respectively. The filled green ( 13 CO ) and yellow ( 12 CO) squares (Oph) and circles (Per) are measurements 
in evenly spaced contour intervals of Ay. The dashed horizontal green and yellow lines are the average of 13 CO and 12 CO contour points. If 
CO traces the mass we find using extinction maps, we would expect the ratio of CO/ Ay to be of order unity as shown by the solid black line. 
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Fig. 6. — S gas and Esfr are show for the sample of massive dense clumps from the survey of IWu et~a l. (2010). Gas surface densities 
are measu r ed from the HCN J=l-0 maps and SFRs are estimated from the t otal I R luminosity, using the extragalactic prescription from 
Kcnnicutt (1998b). The relation between SFR and dense gas from Wu ct al. (2005) is shown (grey solid line) and is extrapolated to lower 
Sgas- We make a cut at Ljr > 10 4,5 Lq, below which the clumps are not massive enough to sample the IMF and lie off a the linear relation 
(Section EH) . 
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Fi g. 7. — Comparison o f massive, dense clumps with Ljr > 10 4 5 Lq to extragalactic relations (Kcnnicutt (1998b), Bigicl et al. (2008), 
and IKrumholz et al.l l|2009D . blue, red, and orange lines, respectively). The relation between SFR and dense gas from lWu et al .1 112003 1 is also 
shown (grey solid line). 
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Fig. 8. — The ratio of Esfr and S gas compared to S gas for low and high-mass star forming regions. We find a steep fall off in £sFR./S ga s 
in the range of E gas ~100-200 Mq pc~ 2 . We denote this steep fall off as a star forming threshold, Ethi between active star forming regions 
and inactive regions. 
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Fig. 9. — Comparison of Galactic total c2d and GB clouds, YSOs, and massive clumps to extragalactic relations. SFR and gas surfaces 
densities for the total c2d and GB clouds (cyan squares), c2d Class I and Flat SED YSOs (green and magenta stars), and Ltr, > lO 4 ' 5 I/0 
massive clumps (y ellow diamonds) are shown. The range of gas surface densities for the spirals and circumnuclear starburst galaxies in the 
IKennicutti jl998bjpl sample are denoted by the grey horizontal lines. The grey shaded region denotes the range for S t h of 129±14 Mq pc . 
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Fig. 10. — We fit Class I and Flat SED YSOs (green stars) and massive clumps (yellow diamonds) to a broken power law (Section 13.211 
and obtain an estimate for the star forming threshold, S t hi of 129±14 Mq pc -2 (grey shaded region). The slope changes from 4.6 below S t h 
to 1.1 above E t h- 



